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FOREWORD 


This  final  technical  report  describes  the  work  accomplished 
during  the  period  24  April  1985  to  30  November  1987  by  the 
University  of  Dayton  Research  Institute  (UDRI)  under  Contract 
Number  DTFA03-85 -C-00009  for  the  Federal  Aviation  Administration. 
Mr.  Ernest  Schlatter  was  the  technical  contract  monitor. 

Cooperation  and  assistance  from  a  number  of  individuals  and 
organizations  was  vitally  important  to  the  success  of  the  program. 
Mr.  P.  Thorson  of  ETAC  and  his  colleagues  provided  contour  map  data 
from  the  ETAC  data  base.  Dr.  R.  K.  Jeck  of  NRL  provided  a  great 
deal  of  input  for  the  collection  of  the  data  as  well  as  providing 
data  for  validating  the  Smith-Feddes  model.  Mr.  C.  W.  Rogers  of 
Calspan  provided  the  computer  code  and  user’s  documents  for  the 
updated  Smith-Feddes  model. 

Ms.  Joanda  D'Antuono  was  responsible  for  typing  and  assembling 
this  report. 
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ABSTRACT 


Realistic  methods  are  discussed  for  identifying  potentialy 
severe  icing  geographical  regions  of  the  Northern  Hemisphere.  A 
useful  method  to  calculate  the  liquid  water  content  of  a  remote 
region  (using  satellite  data)  is  described.  Tne  collection  of  non¬ 
conus  aircraft  is  shown  and  data  is  presented.  ■  ' 
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SECTION  1 
INTRODUCTION 

This  is  the  the  final  report  detailing  the  work  accomplished 
under  Contract  DTFAQ3  -85-C-00009  for  the  Federal  Aviation 
Administration  (FAA).  Essentially  the  contract  requested  research 
in  these  three  areas: 


1^  Survey  foreign  countries  and  retrieve  whenever  possible 
aircraft  icing  data  that  would  be  suitable  for  inclusion 
into  the  already  established  FAA/NRU  icing  data  base. 

2)  Determine  those  geographical  areas  where  icing  would  be 
a  problem  for  aircraft. 

.  ,  '''  3)  Develop  an  alternate  data  collection  scheme  which  could 

be  used  in  lieu  of  flying  aircraft  to  obtain  icing  data. 

1.  1  BACKGROUND  OF  AIRCRAFT  ICING  DATA 

The  current  FAA/NRL  icing  data  base  is  used  to  characterize 
the  potential  icing  environment  thru  which  aircraft  fly.  One  main 
purpose  of  this  data  base  is  to  predict  maximum  icing  conditions, 
so  that  the  proper  cautionary  measures  can  be  taken.  Curently, 
this  data  base  is  primarily  from  the  continental  United  States. 

It  was  the  desire  of  the  FAA  to  broaden  the  data  base  to  include 
data  from  other  parts  of  the  world. 

The  main  problem  in  attaining  icing  data  from  other  coun¬ 
tries  was  that  few  foreign  countries  had  programs  that  utilized 
aircraft  to  collect  data  on  liquid  water  content,  median  volume 
diameter  and  other  parameters  for  the  express  purpose  of  assembly- 
ing  the  data  into  a  data  base  similar  to  that  of  the  FAA's;  France 
and  South  Africa  being  the  only  two. 

Nonetheless  other  countries  such  as  Canada  had  programs  in 
cloud  seeding  from  which  data  could  possibly  be  used  by  the  FAA. 
England  had  three  agencies  which  had  data  that  could  possibly  have 
proven  useful  to  the  FAA,  Australia,  and  Russia  had  possible 
programs,  Germany  had  a  program  that  originally  looked  promising, 
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and  finally  the  JTD  Environmental  Company  located  in  Vasendera, 
California  performed  icing  certification  work  for  foreign  coun¬ 
tries.  These  countries  and  company  were  possible  sources  for  the 
FAA/NRL  icing  data  base.  Data  which  proved  applicable  to  the 
FAA/NRL  data  base  was  obtained  from  the  sources  shown  in  Table 
1.1. 

Four  steps  were  required  for  obtaining  non-continental 
United  States  (non-conus)  icing  data. 

1.  Construct  a  mailing  list  of  all  agencies  thought  to 
have  icing  data.  The  rule  was  -  if  in  doubt  send  a  let¬ 
ter  . 

2.  Use  the  response  of  the  original  letters  for  new  mailing 
or  further  contact  if  the  initial  response  was 

f  avorable . 

3.  When  data  was  located,  determine  the  group  in  the  agency 
and  a  name  in  the  group  who  had  the  authority  to  send 
UDR I  the  data. 

4.  Make  the  proper  arrangements  for  the  data  to  be  sent  to 
UDRI  in  the  desired  format. 

Step  one  was  relatively  easy  due  to  the  work  done  by  Dr. 

Jeck  of  the  NRL  in  alerting  UDRI  to  those  agencies  who  had  data. 
The  results  of  Steps  2,  3,  and  4  are  summarized  by  country  in 
Appendix  A.  Data  that  has  been  received  from  foreign  countries 
that  will  be  considered  for  inclusion  in  the  FAA/NRL  data  base  are 
summarized  in  Table  1.1  of  Section  1.1.  Data  received  in  summary 
form  that  will  not  find  its  way  into  the  data  base  appears  in 
Table  1.2  of  the  same  section.  Appendix  A  list  the  complete  mail¬ 
ing  list  as  well  as  the  responses.  Appendix  B  is  a  replica  of  the 
letters  sent. 

Table  1.1  is  a  listing  in  summary  form,  delineated  by 
country  and  agency,  of  the  data  received  by  UDRI  during  this  con- 
tractural  effort.  Data  was  received  from  Canada  (two  agencies), 
the  JTD  Environmental  Services,  the  South  African  Weather  Bureau, 
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TABLE  1  .  1 


SUMMARY  OF  COLLECTED  FOREIGN  DATA 

DATA  FROM  ALBERTA  RESEARCH  COUNCIL 
ALL  FLIGHTS  MADE  IN  NORTHERN  CANADA 

1983-Winter 


Date 

Temp 

Alt 

MVD 

Max . 

LWC 

Distance/Time 

Type  oft 

Range 

Range 

Range 

J.W. 

FSSP 

in 

clouds 

Clouds 

or' 

kft . 

micron 

grams/m 

km 

min 

m/03 

-10,-14 

09,  1 1 

02,45 

•  37 

.97 

71 

27.3 

SC/CU/TCU 

15/C3 

-13,-15 

1  0,  12 

02,45 

mo 

.84 

96 

35.7 

SC/CU/TCU 

16/03 

-15,-20 

10, 13 

02,45 

•  45 

.71* 

224 

56.8 

SC 

17/03 

-17,-21 

11,13 

02,45 

•  34 

.61* 

63 

20.  3 

SC/AC 

21/03 

-1 1 ,-22 

11,13 

02,45 

.32 

.30* 

25 

15.  3 

AS/ST 

22/03 

-11,-15 

05,11 

04,4  5 

.39 

•  32* 

84 

21  .8 

SC/AS 

23/03 

-01  ,-14 

06,11 

02,45 

1  .60 

1 .  39* 

145 

48.2 

TCU/SC 

24/03 

-09,-18 

07,13 

02,45 

.62 

.93* 

189 

59.2 

CU/TCU 

1983-Summer 


16/05 

-06,-10 

10,11 

02,  31 

.73 

.76* 

98 

28. 3 

CU/TCU 

31/05 

-03,-14 

12,  16 

02,31 

.86 

.66* 

120 

32.0 

CU/TCU 

01/06 

1 

O 

1 

O 

CX> 

00,09 

02,45 

.90 

.80 

4 

6.8 

CU/TCU 

06/06 

O 

I 

=r 

0 

1 

12,15 

02,31 

.89 

.77 

117 

30.2 

CU/TCU 

08/06 

00,-n 

10.15 

02, 31 

1.18 

1.25* 

61 

15.7 

TCU/CB 

09/06 

-10,-12 

14,15 

04,31 

1.00 

1.22 

199 

48.4 

TCU/CB 

10/06 

-04,-13 

12,15 

02,31 

.65 

1.15 

149 

40.0 

CU/TCU 

1  1/06 

-07,-12 

11,13 

02,31 

.74 

1.06 

53 

16.9 

CU/TCU 

12/06 

-01 ,-1  1 

09,13 

02, 31 

.87 

1.18 

62 

18. 3 

CU/TCU 

19/06 

03,-09 

06,09 

02,45 

.46 

.7  3* 

54 

19.2 

ST/AS/NS 

21/06 

-05,-10 

11,13 

05,45 

1.33 

1 .73* 

72 

10. 1 

CU/TCU 

22/06 

-03,-13 

08,  15 

07,45 

1  .17 

2.  03 

48 

12.5 

TCU/CB 

29/06 

-07,-1 1 

14,15 

08,31 

.98 

.90* 

19 

4.8 

CU/TCU 

30/06 

00,-05 

10, 13 

08,31 

.65 

.66 

60 

23. 0 

ST/AS/NS 

03/07 

-07,-1  1 

10, 15 

05,45 

.84 

1.22* 

59 

14.1 

TCU/CB 

07/07 

-09,-11 

16,  17 

02,31 

.  09 

.08* 

1  1 

3.1 

TCU/CB 

08/07 

00,-12 

11,18 

03,  31 

•  53 

.60 

16 

4.5 

TCU/CB 

12/07 

00,-12 

11,16 

02,  31 

1  .50 

1  .70 

109 

25.2 

TCU/CB 

17/07 

-07,-11 

12,15 

02, 31 

1.7  3 

.80* 

75 

18.  1 

TCU/CB 

21/07 

00,-1 1 

10,15 

02,  31 

1  .48 

1  .55 

88 

22.9 

TCU/CB 

26/07 

03.-12 

08, 16 

02,  31 

1.14 

1.27* 

51 

13.  1 

TCU/CB 

28/07 

02,-11 

1  0,  16 

02, 31 

1  .02 

1.12* 

53 

13.3 

TCU/CB 

T  SC-3tratocumuiu3 ,  CU-cumulu3,  AC-altocumulus ,  AS-altoatratus ,  ST-stratus, 
CB-cumulonimbus,  NS-nimbostratus ,  TCU-subset  of  cumulus,  e.g.  cumulus- 
conjectus,  ICON-?  designation  used  by  ARC. 
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Table  1 .  1  (continued) 


Date 

Temp 

Alt 

MVD 

Max. 

LWC 

Distance/Time 

Type  of 

Range 

Range 

Range 

J  .W. 

FSSP 

in 

clouds 

Cl ouds 

°C 

kft . 

micron 

grams/m 

km 

min 

01/08 

-08,-13 

17,  18 

02,  31 

1.29 

1.09 

86 

21  .9 

TCU/CB 

03/08 

-01,-13 

12,18 

02,  31 

1 .  36 

1.34 

77 

21  . 1 

TCU/CB 

14/08 

-09,-13 

15, 16 

02.  31 

.67 

.71* 

51 

1  6.4 

TCU/CB 

15/08 

-06,-12 

15, 16 

02,  31 

.76 

.74* 

88 

18.8 

TCU/CB 

16/08 

-07,-11 

09,16 

02,  31 

1  .22 

1 .07 

66 

15.0 

TCU/CB 

17/08 

-01  ,-1  1 

13,  16 

02,  31 

1.27 

.71* 

99 

21 .4 

TCU/CB 

18/08 

-09,-10 

13,14 

08,  31 

.  07 

.  09* 

26 

5.6 

ST/AS/NS 

24/08 

-02,-10 

12,16 

02,31 

1 .04 

.82 

90 

21.2 

CT/TCU 

28/08 

03,-10 

11,16 

02,  31 

1  .79 

.89 

152 

38. 1 

CT/TCU 

198 

3-Winter 

28/11 

1 

O 

—'J 

1 

(V 

O 

04,11 

02,  31 

.12 

.  19* 

18 

22.7 

No  Information 

08/12 

-15,-20 

11,14 

07,  31 

•  39 

.30 

214 

61.4 

AS/AC 

09/12 

00,-20 

05,07 

01,16 

.97 

.31 

8 

9.0 

AS/AC 

10/12 

-04,-17 

05.13 

01  ,  15 

.98 

.39* 

258 

65.8 

AS/AC 

28/12 

-07,-20 

04,11 

02,31 

.21 

.  18* 

19 

18.0 

No  Information 

1984-Winter 

01/02 

-15,-22 

11,14 

02,31 

.29 

.25* 

138 

34.2 

CU 

02/02 

-21,-22 

14,  1  5 

20,31 

.15 

2.75* 

4 

1.8 

No  Information 

06/02 

-14,-19 

13,15 

12,  31 

.06 

.48* 

7 

1  .9 

AC 

08/02 

-08,-1  0 

11,12 

10,  31 

.37 

.26* 

125 

29.0 

CU 

09/02 

-1 1 ,-22 

11,16 

02,  31 

•  38 

.68 

135 

50.8 

CU 

10/02 

-16,-22 

11,13 

02,  31 

.37 

.44 

259 

83.7 

AC 

11/02 

-03,-15 

04,11 

02,  31 

.57 

.  08 

10 

3.0 

CU 

12/02 

-03,-18 

06, 1  5 

10,  31 

.22 

•  32 

133 

33.9 

CU 

1  3/02 

-09,-19 

11,15 

04, 31 

.  19 

.24 

217 

51.3 

CU 

1984-Summer 

1  3/ 06 

-07,-1 1 

12,14 

02,12 

.99 

.64* 

3^ 

10.  1 

CU/TCU 

15/06 

00,-11 

08,  15 

02,  31 

2.17 

1  .39 

35 

8.6 

TCU/CB 

19/06 

03,-13 

05,13 

02,31 

4.25 

1.52* 

24 

6.7 

CU/TCU 

20/06 

02,-10 

08,  14 

02,  31 

4.02 

1.10* 

60 

15.6 

CU/TCU 

22/06 

01,-11 

07,  14 

02,31 

2.24 

1.10* 

46 

12.0 

CU/TCU 

24/06 

00,-10 

10, 15 

02,  31 

1  .61 

.86 

131 

33.2 

TCU/CB 

25/06 

04,-11 

08,14 

02,  31 

3.  38 

1.31* 

71 

17.7 

TCU/CB 

29/06 

-03,-13 

13,  18 

09,31 

1.16 

1.12* 

17 

5.2 

CU/TCU 

02/07 

-02,-15 

1  0,  is 

02, 31 

1.35 

.62 

38 

10.6 

CU/TCU 

04/07 

04,-12 

11,16 

02,  31 

2.12 

1.10 

105 

24.4 

CU/TCU 

05/07 

00,-12 

1  0, 16 

02,31 

2. 01 

.65 

72 

12.0 

TCU/CB 

09/07 

-01,-11 

10,  14 

02,  31 

1.41 

.79 

99 

25.7 

CU/TCU 

11/07 

-06,-12 

14,16 

02,  31 

1  .98 

1  .80 

147 

29.5 

TCU/CB 

12/07 

05,-1 3 

08, 15 

02,45 

3.58 

1.51* 

178 

37.9 

TCU/CB 
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Table  1 . 1  (continued) 


Date 

Temp 

Alt 

MVD 

Max. 

LWC 

Distance/Time 

Type  of 

Range 

Range 

Range 

J  .w. 

FSSP 

in 

clouds 

Clouds 

°C 

kft. 

micron 

grams/m 

km 

min 

1  3/07 

-06,-10 

12,  14 

02,31 

2.22 

1.15* 

32 

9.9 

CU/TCU 

17/07 

04,-10 

10,  16 

02, 31 

2.80 

1.56* 

106 

21.6 

TCU/CB 

19/07 

-09,-10 

14,16 

12,31 

2.43 

1.10* 

16 

3.4 

TCU/CB 

27/07 

-04,-10 

16,18 

02,  31 

1.14 

.75 

80 

22.7 

CU/TCU 

28/07 

-02,-10 

13.17 

02,31 

3.35 

1 .18* 

174 

43.3 

TCU/CB 

30/07 

01,-01 

12,16 

12,  31 

1.73 

1 .62* 

72 

16.0 

TCU/CB 

03/08 

05,-03 

1  0,  14 

04,08 

.03 

.06 

2 

1.2 

CU/TCU 

04/08 

04,-01 

10, 17 

04,  31 

1.27 

1.60* 

56 

14.8 

CU/TCU 

05/08 

00,-1  1 

13,  17 

02,  31 

1  .95 

1.32* 

65 

20.5 

TCU/CB 

06/08 

00,-12 

13,17 

04,45 

3.42 

1.90 

128 

31.4 

TCU/CB 

09/08 

-05,-12 

15,17 

02,45 

.62 

.61* 

81 

20.5 

CU/TCU 

10/08 

-02,-02 

13.  13 

14,  31 

.31 

.08* 

2 

0.4 

CB 

1  3/08 

-05,-11 

08,15 

02,31 

2.68 

1  .58* 

157 

40. 1 

CU/TCU 

14/08 

04,-05 

07,12 

02,  31 

1.01 

.91* 

68 

20.9 

CU/TCU 

17/08 

-08,-10 

12,  17 

02,31 

1 . 01 

.67 

56 

14.1 

TCU/CB 

18/08 

05,-11 

13,  17 

06,31 

2.69 

1.08 

95 

17.8 

TCU/CB 

20/08 

-05,-15 

11,16 

04,31 

.90 

.52* 

31 

1 1  .0 

CU/TCU 

24/08 

03.-13 

10,16 

02,31 

2.78 

4.77 

1 18 

27.8 

TCU/CB 

25/08 

-05,-10 

12,14 

02,31 

.67 

.59 

88 

17.7 

CU/TCU 

1985-Summer 


03/06 

-02,-10 

09,13 

02,31 

.88 

1 . 18 

41 

12.5 

CU/TCU 

06/06 

-09,-10 

12,13 

02,45 

.57 

.89* 

1  12 

27.4 

CU/TCU 

10/06 

-10,-16 

10,12 

02,45 

.45 

.82* 

42 

11.3 

CU/TCU 

12/06 

-03,-15 

11,17 

02,45 

1.42 

2,04 

176 

43-6 

CU/TCU 

15/06 

-01,-16 

05,  15 

02,31 

.55 

1.35 

178 

50.2 

CU/TCU 

20/06 

00,-11 

08.13 

02, 31 

2.54 

2.57 

137 

35.8 

CU/TCU 

2  3/06 

00,-15 

09,16 

05,31 

2.12 

1.88 

93 

21.0 

TCU/CB 

24/06 

02 , - 1  5 

05,14 

02, 31 

1.39 

1 .69* 

1 18 

32. 0 

CU/TCU 

25/06 

-01,-11 

08,12 

02,31 

.96 

1.28* 

103 

31.0 

CU/TCU 

27/06 

-07,-10 

13, 15 

06,31 

.46 

.60* 

31 

7.8 

CU/TCU 

28/06 

-10,-1 1 

15,16 

08,31 

.79 

1.00* 

23 

6.3 

CU/TCU 

30/06 

-05,-12 

14,  16 

02,31 

.86 

.73* 

38 

10.  1 

CU/TCU 

01/07 

00,-10 

13,16 

02,  31 

2.07 

1.99 

82 

23.0 

CU/TCU 

05/07 

-04,-10 

13, 17 

02,  31 

1.66 

1.68* 

16 

6.2 

TCU/CB 

1  0/07 

00,-12 

14,16 

02,31 

1.25 

1.32* 

1  07 

28.7 

CU/TCU 

1  1/07 

-04,-12 

14,16 

02,31 

2.65 

2.63* 

156 

35.2 

CU/TCU 

12/07 

00,-20 

09,20 

02,31 

2.10 

2.27 

179 

46.6 

CU/TCU 

1  3/07 

03,-13 

09, 14 

02,31 

1.82 

2.41* 

145 

41.0 

CU/TCU 

1  5/ 0 

-05,-12 

12,14 

02,31 

1  .94 

2.64 

108 

35.6 

CU/TCU 

17/07 

03,-12 

06,14 

02,31 

2.54 

2.91* 

110 

31 . 0 

CU/TCU 

18/07 

00,-1 1 

10,15 

02,31 

1  .68 

1  .77* 

113 

27.9 

CU/TCU 

22/07 

01,-16 

11,20 

02,31 

1 .61 

1 .  31* 

155 

43.9 

CU/TCU 

2  3/07 

00,-10 

11,18 

02,  31 

3.25 

3.10 

62 

15.9 

TCU/CB 

24/07 

00, -1  1 

09,16 

02.31 

2.01 

2.34* 

92 

24.0 

CU/TCU 
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Table  1.1  (continued) 


Date 

Temp 

Alt 

MVD 

Max. 

LWC 

Distance/Time 

Type  of 

Range 

Range 

Range 

J.W. 

FSSP 

in 

clouds 

Clouds 

°C 

kft. 

micron 

grams/tn 

km 

min 

29/07 

00,-11 

10,16 

02,08 

1.29 

.40 

209 

61.3 

CU/TCU 

30/07 

“02,-13 

12,16 

01,08 

1.95 

.01 

90 

25. 1 

TCU/CB 

02/08 

02,-11 

11,17 

02,31 

2.06 

2.03* 

137 

34.7 

TCU/CB 

03/08 

-07,-11 

15,  17 

05,  31 

1.90 

1.77* 

1 1 

4 . 1 

TCU/CB 

04/08 

-01,-13 

13, 18 

06,  31 

2.57 

2.42 

96 

22.4 

TCU/CB 

05/08 

-04,-11 

12,14 

02, 3* 

2.42 

2.  34* 

152 

37.3 

CU/TCU 

13/08 

-04,-1  1 

12,14 

05,31 

1.10 

1.23* 

103 

28.3 

CU/TCU 

14/08 

00,-01 

10,  1  0 

10,  31 

.18 

.07* 

9 

2.7 

CU/TCU 

30/08 

-07,-10 

12,15 

02, 31 

2.55 

2.22* 

68 

36.8 

TCU/CB 

DATA  COLLECTED 

BY  THE 

ATMOSPHERIC 

ENVIRONMENTAL  SERVICE,  CANADA 

Flights 

over  London, 

Ontario  in 

1981 

30/09 

01,-07 

05,  10 

08,31 

.19 

.16* 

17 

4 . 0 

SC 

01/10 

01.-13 

01 , 10 

04.31 

1.23 

1.68 

262 

58.4 

sc 

02/10 

01,-07 

02,08 

07/31 

.82 

.52 

321 

72.0 

sc 

15/10 

03,-05 

07,10 

08,31 

.60 

•  53* 

143 

30.9 

SC 

04/1 1 

-05,-1  1 

04,06 

06,31 

.27 

.38* 

20 

4.7 

SC 

05/1  1 

-06,-13 

02,  1  1 

05,31 

1  .61 

2.34* 

166 

39.2 

SC 

06/1  1 

01,-13 

02, 10 

07,31 

1 .72 

.55* 

254 

57.8 

SC 

10/11 

-02,-06 

05,06 

11,31 

.33 

.  18 

85 

18.0 

sc 

16/1  1 

-02,-04 

04,10 

09,  31 

.51 

.25* 

72 

15.5 

sc 

17/1  1 

-05, -06 

01 , 11 

07,  31 

1  .08 

.97* 

247 

56.2 

sc 

19/11 

-04,-05 

03,04 

07,31 

1.07 

.47 

169 

50.0 

sc 

20/1  1 

-03,-09 

00,10 

07,31 

.82 

.96 

389 

94.5 

sc 

Data  collected  at 

North  Bay  Ontario,  1982 

29/06 

-01,-04 

05,07 

10,31 

1  .05 

.66* 

14 

3.5 

SC ,TCU 

30/06 

02,-14 

07,14 

17,31 

2.52 

.03 

29 

7.0 

SC.TCU 

01/07 

-08,-14 

12,15 

09,31 

•  37 

.  ?1* 

13 

3-4 

SC.TCU 

02/07 

05,-12 

02,05 

10,31 

1.48 

.92 

176 

46.5 

SC.TCU 

12/07 

01,-02 

10,  10 

14,31 

.10 

.02* 

4 

1  .0 

sc, rcu 

14/07 

-04,-16 

03,  16 

05,  31 

2.73 

1.61* 

154 

35.0 

SC , TCU 

27/07 

-01,-13 

04,11 

05,  31 

2.97 

1.36 

141 

34.5 

SC.TCU 

Data  from  Syracuse, 

New  York, 

1984 

19/10 

-03,-05 

1  0,  14 

07,  31 

1  .97 

1.27* 

60 

15.2 

CU.ST 

23/10 

00, -01 

05,05 

07,  31 

.23 

.18* 

1 

•  3 

C'J.ST 

29/10 

06,-03 

07,  10 

05,  31 

.89 

.61* 

145 

35.3 

CU.ST 

31/10 

-04,-06 

03,05 

07,31 

.54 

1.76 

232 

64.4 

CU.ST 
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Table  1  . 1  (continued) 


Date 

Temp  Alt 

MVD 

Max. 

LWC 

Distance/Time 

Type  of 

Range  Range 

Range 

J  .W." 

FSSP 

In  clouds 

Clouds 

*C  kft. 

micron 

grams/m 

km 

min 

02/ 11 

-05,-11  03,07 

07.31 

1 .80 

.50 

129 

35.5 

CU.ST 

05/ U 

-10,-13  00,14 

06,31 

.76 

1.06 

296 

84.0 

CU.ST 

06/11 

-07,-13  05,08 

06,31 

.69 

.43» 

96 

25.3 

CU.ST 

09/11 

-03,-11  04,12 

07,31 

.79 

.64 

37  0 

97.8 

CU.ST 

10/11 

-09,-09  03,14 

07.31 

.87 

.66* 

170 

43.7 

CU.ST 

11/11 

-09,-14  02,14 

05,31 

.45 

.43* 

59 

14.4 

CU.ST 

12/11 

-07,-12  01 ,10 

08.31 

.73 

.51* 

178 

48.8 

CU,ST 

13/11 

-04,-09  02,06 

10,31 

.61 

.65 

176 

49.5 

CU,ST 

Flights  over  North  Bay 

Ontario  In 

1984 

20/01 

-20,-23  03,03 

07, 31 

.09 

.05* 

20 

4.5 

CU.ST 

21/01 

-20,-26  03,05 

05.  31 

.01 

.  1  1* 

16 

3.7 

CU.ST 

30/01 

-11  ,-18  03,08 

11 ,31 

.07 

.06 

17* 

4.7 

CU.ST 

01/02 

-13,-16  04,06 

07,31 

.15 

.12* 

8 

2.1 

CU.ST 

02/02 

-08,-09  0  4,06 

09,31 

.15 

.23* 

5 

1.6 

CU.ST 

03/02 

-03,-15  03.09 

06.  31 

.46 

.44 

157 

37.8 

CU.ST 

04/02 

-05,-21  02,  12 

07, 31 

3-09 

.51 

107 

26.3 

CU.ST 

07/02 

-22,-26  09,11 

11,31 

.04 

.08* 

41 

10.1 

CU.ST 

08/02 

-09,-15  03,10 

05,31 

.  12 

.17* 

94 

23.9 

CU.ST 

09/02 

-09,-11  06,08 

17,24 

.17 

.32* 

72 

17.7 

CU.ST 

10/02 

-03,-06  06,08 

05,31 

.22 

.21* 

193 

45.9 

CU.ST 

13/02 

-08,-10  03,14 

05,34 

•  35 

•  37 

77 

17.9 

CU.ST 

14/02 

02,-03  02,05 

06,  31 

.39 

.41 

288 

78.0 

CU.ST 

17/02 

03.-07  03,11 

05,31 

.31 

•39 

284 

71.9 

CU.ST 

19/02 

04,-09  01,11 

05,34 

.28 

.43 

287 

71  .4 

CU.ST 

20/02 

-05,-12  02,03 

07,  31 

.35 

.34 

213 

59.6 

CU.ST 

21/02 

-11,-13  03,04 

07, 31 

.10 

.10* 

125 

34.2 

CU.ST 

24/02 

50,-40  02,19 

07,31 

3.75 

.  15 

1«4 

41.4 

CU ,  ST 

28/02 

-08,-14  02,13 

05,  33 

.09 

,06* 

13 

2.4 

CU.ST 

Flights  over  Halifax,  Nova  Sootla 

in  1986 

28/01 

-04,-04  00,05 

11,31 

.60 

.62 

5 

8.1 

30/01 

-05,-40  00,03 

08,61 

.26 

.29 

1  17 

28.9 

02/02 

-05.-40  00,03 

08,61 

.26 

.29 

15 

5. 1 

05/0? 

00,-40  00,06 

1  1  ,61 

1  .12 

.49 

154 

39.3 

15/02 

-01 ,-08  00,06 

05,61 

.02 

.25 

38 

9.5 

18/02 

-03,-40  00,07 

07,31 

.02 

.17 

274 

71.0 

19/02 

-05,-06  01 ,02 

on,  31 

.02 

.12* 

7 

1.8 

22/02 

00,-12  00,09 

10,33 

.57 

.64 

603 

109.2 

25/02 

-01,-10  01,10 

1  3, 34 

.25 

.49 

41 

10.2 

26/02 

-08,-40  00,09 

05, 33 

.06 

.59 

300 

51.1 

27/0? 

-08,-40  00,05 

04,33 

.12 

.43 

54 

14.3 
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Table  1.1  (continued) 


Date 

Temp 

Alt 

MVD 

Max. 

LWC 

Diatance/Time 

Type  c 

Range 

Range 

Range 

J  .W. 

FSSP 

in 

clouds 

Clouds 

°C 

kft. 

micron 

grains/m 

km 

min 

07/03 

-02,-42 

00,08 

06,34 

.65 

.74 

1 1 1 

56.0 

10/03 

-02,-15 

01,11 

07, 34 

.81 

.26 

261 

59.8 

13/03 

-06,-40 

00, 10 

15,32 

.03 

.  08* 

7 

8.3 

Twin  Otter 

28/01 

-03,-09 

09, 12 

12,33 

.00 

.17* 

12 

2.7 

SC 

29/01 

-07,-15 

02,05 

05,32 

•  37 

.27* 

39 

10.2 

SC 

30/01 

-03.-19 

02,16 

07,42 

.12 

.  18* 

34 

7.4 

SC 

02/02 

-04,-1 1 

01,04 

06,34 

•  22 

.  36* 

12 

3.5 

SC 

05/02 

-02,-13 

01 , 14 

07,31 

1  .72 

.55 

60 

14.4 

sc 

12/02 

-03,-15 

00,05 

03,34 

1.37 

2.69 

653 

170.9 

sc 

14/02 

-07,-15 

00,03 

07,32 

.01 

1  .56 

10 

2.7 

sc 

15/02 

-01  ,-20 

00,  16 

07,  31 

.09 

3,  36 

5 

1  .6 

sc 

18/02 

-02,-05 

04,09 

07,33 

.30 

.40 

156 

37.7 

sc 

22/02 

-04,  -14 

01,14 

08,34 

.29 

.69* 

344 

79.6 

sc 

23/02 

-03.-04 

01,01 

05,  31 

.00 

.02 

1 

0.2 

sc 

24/02 

-04,-15 

02,09 

10,31 

.22 

.49* 

25 

7.  1 

sc 

25/02 

-01  ,-09 

00,05 

07,  33 

.22 

.46* 

29 

8.7 

sc 

26/02 

-20,-22 

12,13 

07,31 

.01 

.08* 

1 

0.3 

sc 

27/02 

-12,-15 

05,06 

07.31 

.  04 

.14 

2 

0.6 

sc 

02/03 

-02,-14 

07,15 

11,33 

.29 

.66 

93 

24.7 

sc 

DATA  FROM 

JTD  ENVIRONMENTAL  SERVICES 

FLIGHTS  MADE  IN  NORTH  SEA 

1984 

-Winter 

22/02 

-06,-07 

02,03 

09,14 

.71 

•  35 

150 

42.0 

ST 

2  3/02 

-09,-10 

03,05 

09,18 

1.21 

.55 

124 

28.0 

ST 

24/02 

-05,-12 

04,  10 

06,28 

.43 

.41 

257 

48.0 

ST 

01/03 

-14,-15 

04,  1  1 

15,30 

.91 

•  36 

197 

30.0 

ST 

19/11 

-06,-13 

05, 10 

08,19 

.52 

.71 

178 

35.5 

ST 

20/  1 1 

-03,-06 

03,06 

13,  19 

.99 

.98 

229 

47.0 

ST 

21/11 

-08,-09 

06,07 

11,20 

.46 

.47 

92 

19.0 

ST 

22/  1 1 

-04,-07 

06,07 

09,  14 

.34 

.32* 

41 

7.5 

CU 

23/11 

-20,-22 

17, 18 

19,  30 

.55 

.41 

121 

20.5 

CU 

06/  12 

-09,-19 

10,  15 

08,21 

.36 

.40 

100 

19.0 

ST 

18/12 

-02,-07 

03,08 

09,  39 

.87 

.78* 

208 

45.7 

ST 

19/12 

-05,-13 

05,08 

12,  17 

.68 

.66* 

160 

34.0 

ST 
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Table  1.1  (continued) 


Date 

Temp 

Alt 

MVD 

Max. 

LWC 

Distance/Time 

Type  of 

Range 

Range 

Range 

J.W. 

FSSP 

in 

clouds 

Clouds 

°C 

kft. 

micron 

grama/m 

km 

min 

1985-Summer 

21/08 

-05,-13 

12,  17 

16,  32 

1.50 

1.41* 

306 

51 .8 

CU 

23/  08 

-04,-08 

10,12 

11,25 

1.79 

1.07 

195 

33*0 

CU 

26/08 

-07,-13 

11,14 

10,27 

1.25 

.88* 

338 

62.0 

CU 

29/08 

-08,-10 

11,13 

11,28 

1.41 

1.10* 

278 

51.0 

CU 

02/09 

00,-12 

10,  14 

12,27 

1.37 

.93* 

93 

18.0 

CU 

04/09 

-09,-1  1 

13,  14 

15,30 

1.09 

1.07* 

31 

5.5 

CU 

1986-W inter 

11/  18 

-05,-07 

08,09 

»7 , 27 

2.22 

.13* 

103 

20.0 

ST/ICON 

FLIGHTS 

MADE  IN  SWEDEN 

1985-Winter 

27/02 

-06,-07 

08,10 

11, 2d 

.47 

.53 

145 

28.0 

ST 

28/ C2 

-07,-1  1 

08, 10 

09,30 

.26 

.22* 

151 

33.0 

ST 

01/03 

-05,-1  1 

03,07 

10,21 

.22 

.41* 

224 

49.0 

ST 

08/03 

-07,-10 

03,04 

10,1  3 

.44 

.47* 

144 

26.0 

ST 

1  3/03 

-03,-05 

02,03 

08, 19 

.49 

.58 

195 

44.0 

ST 

20/03 

-05,-06 

03,04 

07,24 

.69 

.73* 

168 

38.0 

ST 

21/03 

-11,-11 

’0,  10 

10,31 

.40 

.36 

80 

17.  0 

ST 

26/03 

-05.-07 

05,08 

08, 34 

.*•3 

.36 

206 

45. 0 

ST 

05/05 

-04,-13 

05,09 

06,29 

.50 

.37* 

177 

30.0 

CU 

06/05 

00,-06 

05,09 

06,16 

.23 

.  19 

83 

12.5 

ST 

1  2/05 

-08,-19 

13,17 

07,27 

.35 

.66 

246 

45.0 

CU 

13/05 

-05,-12 

13,  17 

17,28 

.10 

.04* 

218 

41.0 

ST 

FLIGHTS 

MADE  IN  : 

SPAIN 

1986-Winter 

30/01 

-14,-19 

09,11 

11,28 

.48 

.40 

59 

10.0 

CU 

01/02 

-15,-16 

12,13 

18,27 

.87 

N/O 

77 

13.0 

CU 

12/02 

-07,-08 

11,12 

08, 14 

.07 

.02* 

16 

2.5 

CU 

FLIGHTS 

MADE  IN  IRELAND 

1986-Winter 

1  1/02 

-05,-06 

04,05 

N/0 

.50 

N/O 

198 

39.0 

ST 
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Table  1.1  (continued) 


Oate 

Temp 

Alt 

MVD 

Max. 

LWC 

Distance/Time 

Type  of 

Range 

Range 

Range 

J.W. 

FSSP 

in 

clouds 

Clouds 

°C 

kft . 

micron 

grama/ra 

km 

min 

FLIGHTS 

MADE  IN 

HOLLAND 

1986-Winter 

17/11 

-05,-05 

08,08 

14,20 

.22 

.01 

5 

1 . 0 

ST 

FLIGHTS 

MADE  IN 

NORWAY 

1986 

26/11 

-10,-15 

08,11 

25,33 

1  .66 

.64 

135 

26.0 

CU 

28/11 

-05, -07 

09,11 

05,  30 

.24 

•  31 

153 

32.0 

ST/ ICON 

29/11 

-10,-13 

07.08 

17,30 

1 . 04 

.50 

180 

36.0 

CU 

DATA  FROM 

SOUTH 

AFRICAN 

WEATHER 

BUREAU 

FLIGHTS  OVER  SOUTH  AFRICA 
1984-Winter 


02/10 

-10,-11 

05,05 

07.09 

1.00 

.01* 

12 

2.  1 

CU 

03/10 

-09,-10 

05,05 

07,11 

1.50 

.01* 

21 

3-5 

CU 

11/10 

-12,-14 

05,05 

11,21 

2.70 

.05 

32 

5. 1 

CU 

15/10 

-14,-17 

C5.05 

10,11 

1 .  30 

.10* 

7 

1 . 1 

CU 

17/10 

-04,-14 

05,05 

07,08 

.90 

.10* 

4 

.7 

CU 

18/10 

-10,-12 

05,05 

07, 15 

1.40 

.30 

28 

4.4 

CU 

19/10 

-09,-13 

05.05 

09,10 

.07 

.10 

13 

2.  1 

CU 

22/10 

-12,-15 

05,05 

08,22 

1.50 

.50 

30 

4.7 

CU 

2  3/10 

-12,-13 

05,05 

09,15 

.80 

.20 

20 

3.0 

CU 

24/10 

-09,-12 

05.05 

04,13 

2.  10 

.50 

29 

4.1 

CU 

30/10 

-10,-10 

05,05 

06,06 

.01 

.00 

0 

0 

CU 

07/11 

-13,-14 

05,05 

10,1  1 

1  .00 

.00 

2 

.3 

CU 

12/11 

-11,-12 

05,05 

09,10 

.60 

.00 

1 

.  1 

CU 

13/11 

-03,-13 

04,05 

09, 13 

1 .40 

.00 

1 

.2 

CU 

14/11 

-13,-13 

05,05 

12,12 

2.10 

.00 

1 

.  1 

CU 

26/11 

-09,-17 

05,06 

09,12 

1  .80 

1 .  30* 

14 

2.4 

CU 

27/11 

-10,-15 

04,05 

08, 13 

3.  00 

1.30* 

48 

6. 1 

CU 

07/12 

-09,-14 

05,05 

07, 1  1 

1.10 

.50* 

9 

1.3 

CU 

11/12 

-12,-17 

05,06 

10,  12 

1  .50 

.80* 

8 

1.5 

CU 

12/12 

-04,-08 

05,05 

11,15 

2.  50 

1.10* 

5 

.7 

CU 

13/12 

-11,-25 

05,08 

07, 19 

3.  00 

.90 

20 

3.4 

CU 

14/12 

-07,-19 

04,07 

07,27 

2.40 

.89 

8 

1.3 

CU 

17/12 

-12,-16 

06,07 

09,14 

1  .90 

.90 

16 

2.5 

CU 

18/12 

-10,-17 

05,06 

07,14 

2.30 

1.30* 

30 

4.9 

CU 

19/12 

12,-15 

06 , 06 

06,15 

2.  30 

1 .50* 

26 

3.9 

CU 

20/12 

-12,-18 

06,07 

12,  16 

2.50 

1 . 30* 

14 

2.3 

CU 

21/12 

-08,-23 

06,08 

09,26 

2.70 

.70 

18 

2.8 

CU 

Table  1.1  (continued) 


Date 

Temp 

Alt 

MVD 

Max. 

LWC 

Distance/Time 

Type  of 

Range 

Range 

Range 

J.W. 

FSSP 

in 

clouds 

Clouds 

°C 

kft. 

micron 

grams/m 

km 

min 

1985-WINTER 

02/01 

-03.“08 

04,05 

04,24 

2.50 

.40 

33 

5.1 

CU 

03/01 

-14,-19 

06,07 

08,21 

2.50 

.50 

19 

3.0 

CU 

OH/ 01 

-10,-16 

06,06 

07,30 

1.20 

•  30 

10 

1.5 

CU 

07/01 

-10,-16 

05,05 

07,1  1 

1 .60 

.4  0* 

15 

2.2 

CU 

10/01 

-10,-19 

05,06 

05,22 

2.60 

.50* 

32 

4.9 

CU 

11/01 

-1 1,-16 

05,06 

07, 10 

2.30 

.06* 

18 

3.0 

CU 

17/01 

02,-10 

04,05 

07,  18 

3.10 

.90* 

16 

2.5 

CU 

18/01 

-12,-15 

07,07 

10,26 

3. 1  0 

1.50* 

31 

5.  1 

CU 

21/01 

-01,-14 

06,07 

08, 12 

2.90 

.70* 

18 

2.6 

CU 

22/01 

-10,-12 

06,06 

09,12 

3.00 

.80* 

27 

4.0 

CU 

24/01 

-09,-12 

06,06 

09,  12 

2.50 

.7  0* 

1 1 

1  .7 

CU 

28/01 

-10,-12 

06,06 

10,25 

2.20 

.60* 

21 

3.0 

CU 

29/01 

-10,-15 

06,07 

07,12 

2.70 

.90 

27 

4.2 

CU 

30/01 

-1 1,-14 

06,07 

03,07 

5.00 

.  10 

22 

3.5 

CU 

05/02 

-10,-19 

07,08 

04,08 

2.20 

.20* 

29 

3-3 

CU 

11/02 

-01 ,-08 

05,06 

10,  15 

2.80 

2.30* 

13 

2.0 

CU 

1  3/02 

-02,-11 

04,05 

04,  18 

2.30 

2.10* 

33 

5.4 

CU 

14/02 

08,-05 

03,04 

08,13 

1.20 

1.30* 

9 

1.5 

CU 

15/02 

-12,-18 

05,06 

06,21 

3.00 

2.30* 

82 

8.5 

CU 

18/02 

-09,-16 

06,07 

06, 16 

2.60 

1 .70* 

42 

6.1 

CU 

21/02 

-10,-21 

06,07 

04,15 

2.50 

2.30 

33 

5.3 

CU 

22/02 

-11 ,-14 

06,06 

09,1  3 

2.30 

1.50* 

10 

1  .4 

CU 

25/02 

-09,-12 

06,06 

08,13 

2.4  0 

1 .50* 

30 

4.8 

r.u 

26/02 

-1  1,-15 

06,07 

07, 14 

2.70 

2.50* 

19 

3.1 

CU 

27/02 

-1  0,-14 

06,06 

09,  16 

2.70 

2.00* 

59 

9.0 

CU 

28/02 

-01,-14 

04,06 

06,15 

2.30 

2.6C* 

21 

3-5 

CU 

01/03 

-09,-10 

05,05 

11,16 

2.50 

1  .90* 

5 

.8 

CU 

07/03 

-12,-13 

05,05 

06,13 

2.10 

1.30* 

41 

6.7 

CU 

08/03 

-1  1,-12 

05,05 

08,13 

1  .80 

1.30 

12 

1  .9 

CU 

11/03 

-11,-13 

05,05 

11,15 

2.50 

1 .70 

6 

1.0 

CU 

12/03 

-08,-11 

05,05 

1M3 

2.20 

1.50* 

18 

2.8 

CU 

13/03 

01,-04 

04,05 

05,1  3 

1  . 50 

1.30* 

22 

3.8 

CU 

18/03 

-12,-15 

05,05 

05, 15 

2.90 

2.20 

13 

2.  1 

CU 

26/03 

-03,-16 

04,05 

08, 15 

2.00 

1.50* 

16 

2.4 

CU 

DATA  COLLECTED  BY  THE  LABCRATOIRE  ASSOCIE  DE  METEOROLOGIE  PHYSIQUE 
All  Flights  over  Southern  France-1984 


05/06 

-04,-21 

02,08 

07,35 

2.58 

4.59* 

82 

15.  0 

CU 

06/06 

-04,-21 

02,08 

08,35 

1 .89 

1.80 

147 

26.7 

CU 

07/06 

-00,-27 

02,08 

08,  36 

1 .95 

2.33* 

74 

12.9 

CU 

1-11 


Table  1.1  (continued) 


Date 

Temp 

Alt 

MVD 

Max. 

LWC 

Distance/Time 

Type  of 

Range 

Range 

Range 

J.W. 

FSSP 

in 

clouds 

Clouds 

°C 

kft. 

micron 

grams/m 

km 

min 

19/06 

-06,-21 

02,07 

07,36 

0.00 

1.07 

1  19 

20.6 

CU 

26/06 

-07,-13 

02,04 

07,3*1 

1.21 

1 .61 

169 

26.8 

CU 

27/06 

-00,-24 

05,07 

07,38 

1  .44 

1.09* 

23 

3.8 

CU 

28/06 

-05,-18 

04,08 

08,35 

1  .70 

2.76 

115 

20.0 

CU 

01/07 

-09,-23 

06,08 

08,32 

1 .40 

1 .64 

85 

1  3-2 

AS 

02/07 

-05,-07 

04,06 

05,16 

.76 

.44* 

20 

3.2 

ST 

DATA  FROM  THE  UNIVERSITY  OF  CAMBRIDGE,  ENGLAND 


This  data  is  different  than  flight  data  in  that  it  is  data  from  a  tower  which 
is  enveloped  in  clouds  for  long  periods  of  time.  Thus  there  is  no  range  in 
altitude  or  distance  of  penetration.  The  following  is  five  days  of  cloud 
data  and  two  days  of  ground  fog  data. 


06/04 

07/04 

12/04 

01/-02 

00/-01 

01/-02 

06,24 
05,21 
13,  16 

Cloud  Data 
1984 

1.**3 

.96 

1.33 

97 

115 

1  30 

1985 

1  1/08 
12/12 

01/-02 
00/ -01 

16,  17 
11/12 

.60 

.51 

125 

1 1 

1982 

Ground  Fog 

Data 

24/03 

00/-01 

10/12 

.57 

108 

1983 

Ground  Fog 

Data 

27/10 

00/-01 

10/12 

.85 

123 
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TABLE  1.2 


SUMMARY  OF  GERMAN  DATA 


GERMANY 
1983/  1984 


Date 

Alt. 

km 

Temp. 

°C 

LWC 

gm/m 

9/1 

02 

-14 

.  18 

16/  1 

03 

-12 

.  17 

17/  1 

02 

-03 

.36 

18/  1 

02 

-13 

.37 

25/1 

02 

-12 

.47 

2/3 

02 

-06 

.14 

8/3 

02 

-12 

.09 

9/3 

02 

-13 

.25 

27/3 

03 

-10 

.  19 

28/3 

03 

-06 

.  18 

29/3 

03 

-09 

.17 

3/4 

02 

-1 1 

.37 

10/12 

03 

-03 

.09 

11/12 

03 

-03 

.10 

25/01 

02 

-10 

.24 

30/01 

03 

-05 

.13 

06/02 

02 

-04 

.36 

07/02 

02 

-05 

.79 

1  1/03 

02 

-1 1 

.32 

12/03 

01 

-07 

.43 

26/04 

01 

-03 

.30 

29/04 

02 

-10 

.44 

05/03 

03 

-04 

.08 

07/03 

02 

-05 

.21 

13/03 

01 

-04 

.49 

21/03 

02 

-1  1 

.15 

24/03 

03 

-03 

.19 

03/  04 

03 

-06 

.07 

04/04 

02 

-05 

.30 

MVD 

micron 

Dtst. 

km 

Time 

min 

Cloud 

- 

297 

75 

ST,  SC 

23 

233 

58 

CU 

27 

600 

149 

AS,  CU, 

21 

440 

108 

SC,  ST 

15 

559 

137 

AC,  AS 

13 

197 

49 

AC,  AS 

- 

201 

50 

AC,  AS 

- 

318 

78 

AC,  AS 

15 

583 

145 

AC ,  AS 

25 

257 

64 

AS 

- 

297 

74 

AS 

17 

688 

170 

SC,  AC 

1984/1985 

33 

257 

63 

ST 

25 

359 

88 

SC 

19 

478 

1  19 

AC 

36 

352 

89 

AC,  AS 

- 

276 

71 

AC,  AS 

21 

663 

170 

NS,  AC 

15 

633 

165 

SC,  AC 

21 

4  05 

104 

ST,  AS 

- 

299 

79 

ST,  AS 

15 

363 

93 

AC 

1985/  1986 

27 

397 

92 

ST 

27 

487 

123 

ST 

23 

290 

78 

ST 

19 

367 

90 

AS,  AC 

21 

290 

70 

AC 

13 

493 

1 1 1 

AS,  AC 

17 

322 

214 

SC 

NS 
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France  and  England.  The  data  received  was  normally  collected  at 
one  hertz  intervals  and  included  altitude,  ambient  temperature, 
liquid  water  content  from  two  instruments,  median  volume  diameter, 
and  cloud  type. 

In  some  cases  (particularly  the  Canadian),  the  incoming  data 
included  much  time  during  which  the  aircraft  was  not  flying  in  a 
cloud.  Hence  the  incoming  data  was  filtered  through  a  computer 
program  that  eliminated  these  data  and  also  averaged  the  one  hertz 
data  over  one  minute.  This  data  was  then  subjected  to  visual 
scrutiny  for  obvious  instrument  error  or  other  types  of  malfunc¬ 
tions  . 

Finally,  this  averaged  data  was  further  summarized  by  deter¬ 
mining  given  values  for  each  day  of  flight.  These  are  the  values 
that  appear  in  Table  1.1: 

1  Date 

2  Temperature  range,  °C 

3  Altitude  range,  kft 

4  Median  volume  diameter  range,  micrometers 

5  Maximum  Liquid  Water  Content  from  both  the  LWC 
and  FSSP.  The  asterik  signifies  that  the  maxi¬ 
mum  reading  from  both  instruments  occurred  at 
the  same  time. 

6  Kilometer  and  minutes  in  clouds 

7  Type  of  clouds. 

Table  1.1  is  intended  to  give  the  reader  a  sense  of  the 
quality  and  quantity  of  data  received. 

1.2  NON-MAGNETIC  TAPE  DATA 

Unfortunately,  the  data  received  from  Germany  was  not  one 
hertz  data.  Also,  it  was  received  on  hardcopy  printout  in  tabular 
form  and  all  attempts  to  acquire  the  original  aircraft  data 
failed.  However,  a  summary  of  the  data  received  from  Germany  is 
shown  in  Table  1 .  2. 
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SECTION  2 

GEOGRAPHICAL  OCCURRENCE  OF  SEVERE  ICING 


This  section  discusses  the  world-wide  occurrence  of  "severe" 
icing.  Later  in  the  section  a  workable  definition  of  severe  icing 
will  t>e  developed  in  terms  of  characterizing  the  world-wide  icing 
environment.  The  occurrence  of  severe  icing  is  much  more 
widespread  than  is  shown  by  the  existing  or  future  FAA  icing  data 
base.  Therefore,  there  is  no  guarantee  that  a  characterization  of 
the  icing  atmosphere  based  on  the  available  data  is  either  repre¬ 
sentative  or  conservative.  A  first  step  to  remedy  this 
shortcoming  is  to  identify  the  extent  of  severe  icing. 
Unfortunately  this  is  not  currently  feasible  on  a  global  scale, 
however,  it  is  possible  for  a  large  part  of  the  northern  hemi¬ 
sphere  primarily  north  of  20  to  30°  north  latitude.  After 
identification,  the  areas  of  severe  icing  for  which  there  are 
measurements  can  be  compared  to  those  without  measurements.  From 
this  comparison  two  conclusions  can  be  drawn.  One,  the  extent  of 
severe  icing  without  measurements  is  much  greater  than  that  with 
measurements.  Two,  it  is  likely  that  the  most  severe  icing  en¬ 
vironment  lies  within  the  former  areas.  These  conclusions  will  be 
expanded  upon  with  more  evidential  support  later. 

2.1  DEFINITION  OF  "SEVERE"  ICING 

It  is  necessary  to  define  "severe"  icing.  A  few  philosophi¬ 
cal  comments  are  apropo.  The  definition  has  to  be  meaningful  in 
aviation  terms  ( 1 . e . ,  one  would  expect  most  aircraft  to  have  dif¬ 
ficulty  in  the  event  of  a  "severe"  icing  encounter);  however,  the 
definition  has  to  be  simple  enough  to  be  applicable  and  has  to 
reflect  existing  technology  with  regard  to  icing  parameters  (see 
Section  3.1).  One  difficulty  is  that  aviation  and  meteorology  use 
different  parameters  and  units  for  describing  icing,  many  of  which 
are  infrequently  measured.  But  the  core  of  the  problem  is  that 
the  capability  to  predict  the  amount  and  nature  of  ice  accretion 
given  the  simple  problem  of  constant  ambient  conditions  has  not 
yet  been  acceptably  established.  And  finally  the  realistic 
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problem  of  variable  cloud  conditions  has  not  been  satisfactorily 
addressed.  Consider  previous  definitions  of  severe  icing. 

The  first  consideration  concerns  temperature.  Obviously, 
temperature  must  be  below  0°C  for  freezing  to  occur.  Depending  on 
the  type  of  aircraft  and  its  operation,  however,  the  outside  tem¬ 
perature  must  be  several  degrees  below  freezing  for  icing  to 
occur.  Nonetheless  it  is  not  feasible  to  account  for  all  the  pos¬ 
sible  combinations  of  aircraft  and  usage  so  we  shall  simply  adhere 
to  the  assumption  that  the  temperature  must  be  below  0°C  for  ice 
to  form.  Additionally,  we  assume  that  no  ice  will  form  for  tem¬ 
perature  below  - 4 0 ° C  as  no  liquid  water  is  expected  at  this 
temperatu..  e  and  lower  temperatures. 

The  second  consideration  involves  the  amount  of  cloud  water 

or  liquid  water  content  (LWC).  Generally,  other  things  being 

equal,  a  higher  LWC  implies  greater  ice  accretion.  The  nature  of 

the  ice  accretion  depends  on  the  temperature,  but  as  previously 

mentioned  is  not  yet  acceptably  predictable  for  all  temperatures. 

Therefore,  one  cannot  now  formulate  a  severe  icing  definition  that 

considers  LWC  as  a  function  of  temperature  so  that,  for  example,  a 
3 

LWC  of  .5  gm/m  at  -5°C  would  be  considered  as  severe  as  a  LWC  of 
1  gm/m^  at  a  temperature  of  -15°C.  In  the  future,  such  a  defini¬ 
tion  would  be  desirable  and  would  certainly  make  the  definition 
more  realistic  in  terms  of  the  severity  an  aircraft  would  actually 
experience.  The  popular  LWC  value  for  severe  icing  seems  to  be  1 
gram  of  LWC  per  1  meter  volume  of  air  with  no  consideration  of 
temperature. 

The  medium  volume  diameter  (MVD)  drop  size  is  another  factor 
that  influences  the  severity  of  icing.  Larger  drops  having 
greater  inertia  are  not  only  more  likely  to  strike  an  aircraft  in¬ 
creasing  the  resulting  ice  accretion  but  also  strike  the  aircraft 
on  unprotected  areas.  Thus,  as  far  as  aircraft  are  concerned,  the 
MVD  is  very  important  to  severe  icing.  Unfortunately,  it  is  not 
possible  to  delineate  the  global  icing  environment  in  terms  of  the 
MVD  because  as  discussed  in  Section  3.1  there  are  a  plethora  of 
environmental  factors  which  influence  the  drop  size  distribution 
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and  hence  the  MVD .  In  order  to  specify  the  MV D  for  given  environ¬ 
mental  conditions,  it  is  at  least  necessary  to  have  a  priori 
knowledge  of  these  environmental  factors.  Because  it  is  currently 
impossible  to  do  this  we  are  forced  to  ignore  the  important  MVD 
value  in  our  derivation  of  a  "severe"  icing  definition. 

A  final  word  on  horizontal  and  vertical  extent  of  severe 

•3 

icing  is  also  in  order.  Obviously,  moderate  LWC's  below  1  gm/m 
extending  horizontally  for  great  distances  and  of  considerable 
vertical  extent  are  more  serious  to  an  aircraft  than  are  spatially 
limited  occurrences  of  higher  LWC's.  Isolated  instances  of  ex¬ 
tremely  high  LWC  such  as  occurring  in  cumulonimbus  clouds  are 
rej,J  'y  avoided  by  an  aircraft  while  horizontally  extensive  layers 
of  ...^derate  LWC  are  not  easily  avoided.  This  is  why  a  severe 
icing  definition  of  LWC  >  1  gm/ra^  for  -40°C  <  T  <  0°C  is  sensible 
since  it  applies  to  the  extreme  values  of  layer  cloud  liquid 
water . 

2.2  SEVERE  ICING  ENVIRONMENT 

2.  2.  1  Antecedent  Studies 

The  first  task  in  characterizing  the  global  severe 
icing  environment  is  to  identify  regions  where  severe  icing  occurs 
frequently.  This  task  though  is  complicated  by  a  lack  of  directly 
measured  data  for  regions  where  the  most  frequent  severe  icing 
conditions  are  expected  (i.e.,  oceanic  areas  to  the  east  of 
continents).  Therefore,  it  is  necessary  to  resort  to  indirect 
methods  such  as  numerical  calculation  of  LWC.  Numerical  model 
results  can  be  corroborated  where  direct  or  indirect  measurements 
are  available.  Before  considering  numerical  model  results  it  is 
useful  to  briefly  review  earlier  work. 

During  the  1950’s,  the  Air  Weather  Service  flew  a 
number  of  icing  reconnaissance  flights  over  the  North  Pacific  and 
North  Atlantic  Oceans.  The  reconnaissance  aircraft  flew  at  the 
same  levels  (700  mb  and  500  mb)  on  the  same  routes  everyday. 
Therefore,  it  was  possible  to  statistically  analyze  the  results 
(Appleman,  1959)  in  the  form  of  a  conditional  probability  of 
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aircraft  icing  given  cloud  amounts  >  6/10  as  a  function  of  tem¬ 
perature  and  altitude. 

Katz  (1967)  used  these  results  in  addition  to  tem¬ 
perature  and  cloud  data  to  calculate  the  probability  of  icing  in 
the  intervals  C  to  5,000  feet,  5,000  to  10,000  feet,  10,000  to 
15,000  feet,  and  15,000  to  20,000  feet  by  season  over  the  northern 
hemisphere.  He  assumed  that  the  mean  cloud  amount  at  an  altitude 
is  equal  to  the  probability  of  cloud  amounts  >  6/10  at  that  same 
altitude.  As  a  result  his  probability  of  icing  is 

P ( I )  -  ;“2°  F[ I ( T )  /  C  >  6/10]  P[C  >  6/10]  F(T)  dt 
-40  °C 


where 


F C I ( T ) J  '  0  >  6/10]  -  conditional  probability  density  func¬ 
tion  of  icing  given  a  cloud  amount 
6/10  as  a  f'  notion  of  temperature, 

P[C  >  6/10]  -  probability  of  cloud  amount  greater  or 
equal  to  6 /  10 ,  and 

F(T)  -  probability  density  function  of  temperature. 

The  results  of  Katz’s  study  for  0  to  5,000  feet  in 
winter  give  the  highest  probabilities  of  icing  occurrence  over 
eastern  Canada  and  southern  Greenland  P  >  0. 15,  northern  central 
Europe  P  >  0.15,  and  the  western  and  northern  Pacific  p  >  0.10. 

In  spring  the  dist  ibution  is  similar  except  for  P  >  0. 15  over  the 
northern  Pacific.  In  summer  and  fall  the  probabilities  are 
generally  less  except  over  the  Bering  Strait  i.e.,  P  <0.15  in  the 
fall.  In  general,  the  probabilities  are  less  with  increasing  al¬ 
titude  with  the  exception  of  southeast  Asia  during  the  summer 
monsoon  months,  at  15,000  to  20,000  feet  P  >  0.10.  For  our  pur¬ 
poses,  the  Katz  study  has  two  major  deficiencies.  One,  the  study 
is  based  on  limited  data  and  therefore  may  not  be  truly  repre¬ 
sentative  for  those  areas  without  much  data  (i.e.,  no 
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rawlnsondes ) ,  and  two,  Katz  doesn't  really  consider  "severe" 
icing . 

In  Heath  (197  2),  the  frequency  of  icing  occurrence 
was  calculated  for  the  northern  hemisphere.  Temperature  and  dew¬ 
point  data  for  380  radiosonde  stations  were  analyzed  for  frequency 
of  occurrence  of  temperature/dew-point  temperature  differences 
when  the  temperature  is  below  0°C  for  the  1,000  mb,  350  mb,  700 
mb,  and  500  mb  pressure  surfaces.  These  differences  were  corre¬ 
lated  to  a  probability  of  icing  based  on  the  results  from  the  Air 
Weather  Service  reconnaissance  flights.  The  final  icing  probabil¬ 
ity  was  then  an  integral  over  temperature  of  the  probability 
density  function  of  temperature/dew-point  temperature  difference 
times  the  probability  of  icing  given  those  conditions.  Heath's 
results  show  that  a  high  frequency  of  icing  occurrence  is  expected 
over  southeastern  Canada  at  850  mb  P  >  0.10,  over  Scandinavia  and 
northeastern  Russia  P  >  0.15,  and  over  the  northern  and  western 
Pacific  P  >  0.20  in  January.  Again,  as  for  the  Katz  results, 
there  is  not  enough  data  to  guarantee  full  resolution  of  icing  oc¬ 
currence,  and  again  severe  icing  is  not  considered.  In  addition, 
Katz  and  Heath  somewhat  differ  in  their  results  although  the  over¬ 
all  fields  are  similar. 

Other  studies  have  been  done  but  are  generally 
regional  or  local  in  nature.  Roach  et  al.  (  198  4)  analyzed  20,000 
World  War  II  reconnaissance  flights  along  with  some  recent  flights 
in  producing  a  supercooled  cloud  climatology  for  the  northeastern 
Atlantic  Ocean.  The  analysis  shows  that  the  probability  of  super¬ 
cooled  cloud  with  LWC  i  0.5  gm/ra^  exceeds  0.02  at  850  mb  during 
winter  for  an  area  to  the  northwest  of  the  British  Isles  and  south 
of  Iceland.  Probabilities  for  exceeding  this  threshold  are  lower 
elsewhere  and  at  all  other  vertical  levels  analyzed.  The  prob¬ 
ability  for  exceeding  0.2  gm/m^  at  850  mb  during  winter  for  thp 
same  area  is  0.10.  The  latter  figure  is  consistent  with  the  prob¬ 
ability  of  icing  found  by  Heath  and  similar  to  the  value  found  by 
Katz.  We  know  of  no  similar  studies  for  the  southern  hemisphere. 
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2.3  ETAC  DATA  BASE  MAPS 


In  1984,  the  United  States  Air  Force  Environmental  Technical 
Applications  Center  (ETAC)  compiled  an  extensive  liquid  water  con¬ 
tent  (LWC)  data  base  for  the  years  1977-1980.  In  1986  the  LWC 
values  were  calculated  on  the  Air  Weather  Service  (AWS)  northern 
hemisphere  IJ  grid  for  much  of  the  northern  hemisphere  (north  of 
30°  north)  and  at  15  vertical  levels.  The  approximate  horizontal 
spacing  of  these  grid  points  is  200  nautical  miles;  somewhat  less 
for  lower  latitude  grid  points  and  somewhat  more  for  higher 
latitude  points.  The  original  Smith-Feddes  model  (Smith  1 97  H . 
Feddes  1974)  described  in  Section  3.2  was  used  with  archived  twice 
daily  three-dimensional  nephanalyals  (3DNEPH)  and  AWS  analyses  for 
the  above  mentioned  4-year  period.  3DNEPH  analyses  were  produced 
at  a  horizontal  resolution  of  25  nautical  miles  so  that  an  LWC 
grid  point  is  representative  of  64  3DNEPH  values.  As  a  result, 
approximately  240  LWC  values  were  generated  at  each  vertical  level 
at  points  north  of  30°  on  the  IJ  grid  for  eaoh  month  of  the  year. 

3 DNEPH  analyses  of  tropical  and  southern  hemisphere  regions  were 
apparently  lacking  therefore  the  LWC  values  for  theoo  areas  could 
not  be  calculated. 

The  conditional  probability  of  LWC  values  greater  than  1 
gm/m  with  temperatures  between  0°C  and  -40°C  was  calculated  an¬ 
nually  and  monthly  from  the  data  base.  This  probability  was 
calculated  in  the  eight  vertical  layers  listed  below  in  Table  2.1 
at  each  IJ  grid  point. 

TABLE  2.  1 

SEVERE  ICING  PROBABILITY  CALCULATION  LEVELS 


Layer 

Altitude  Range  (m) 

1 

1  - 

500 

2 

50  1  - 

1000 

3 

1001  - 

1500 

4 

150  1  - 

2000 

5 

200  1  - 

2500 

6 

250  1  - 

3000 

7 

3001  - 

4500 

8 

4501  - 

6000 
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The  oondltlonal  probabilities  are  indicative  of  the  fre¬ 
quency  of  ocourrenct  of  severe  icing  as  defined  in  Section  2.1 
However,  the  Sol th-Feddes  model  probably  exhibits  a  systematic 
bias  for  overestimating  LWC  values.  Therefore,  while  the  result- 
lng  probacil ities  of  LWC  greater  than  1  gm/mJ  are  high,  the 
distribution  of  higher  probabilities  properly  delineates  the 
severe  icing  environment. 

The  monthly  variation  of  the  results  are  similar  within 
seasons.  For  example,  the  results  for  December  and  February  are 
quite  similar  to  those  for  January  arid  results  for  June  and  August 
are  quite  similar  to  those  fer  July,  etc.  Thus,  the  months  of 
January,  April,  July,  and  October  are  representative  of  the 
winter,  spring,  summer,  and  fall  seasons,  respectively. 

F  igut  ••  2.1  shows  the  distribution  of  '’severe  icing"  for  the 
seoond  layer  (501-1000  m)  during  April.  In  general,  probabilities 
for  layer  1  (  1-500  m)  are  quite  low  as  the  cloud  base  is  fre¬ 
quently  higher  than  500  m  and  higher  LWC  valuos  usually  occur  in 
proximity  to  the  cloud  top  rather  than  cloud  base.  Higher  prob¬ 
abilities  OIOS)  are  limited  to  higher  latitudes  with  the  most 
extensive  occurrences  near  and  south  of  the  Kamchatka  peninsula, 
over  eastern  Siberia  from  Lake  Baykal  to  Sakhalin  island,  north  of 
the  Yenlsey  in  the  northern  Soviet  Union,  over  the  north  Atlantic 
from  Greenland  to  the  North  Sea,  over  Norway  and  Sweden,  over  part 
of  central  Europe,  over  southern  Alaska,  and  over  eastern  Quebec 
and  Labrador.  Especially  high  probabilities  are  located  just 
north  and  east  of  Norway  and  over  the  Kunlun  mountain  range  north 
of  Tibet.  This  latter  feature  appears  in  many  of  the  icing  prob¬ 
ability  charts  and  may  be  spurious. 

Figure  2,2  shows  the  distribution  for  layer  4  (150  1  -  2000 
m)  in  April.  In  general,  the  probabilities  are  higher  and  the 
severe  icing  occurrence  is  more  extensive  than  for  layer  2.  The 
oaourrenoe  is  particularly  more  widespread  from  the  Great  Lakes 
northeastward  through  northern  New  England  and  Quebec.  There  is 
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Figure  2.1  Occurrence  of  Supercooled  Liquid  Water  Content  Greater 
than  Igm/m1  in  April  at  501-1000  meters.  (Map  shows 
only  those  regions  where  the  frequency  of  occurrence  is 
10$  or  greater . ) 
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Occurrence  of  Supercooled  Liquid  Water  Content  Greater 
than  Igm/m1  in  April  at  1501-2000  meters.  (Map  shows 
only  those  regions  where  the  frequency  of  occurrence  is 
10$  or  greater . ) 


F igure  2 . 2 


also  another  prominent  area  located  over  British  Columbia  that 
does  not  appear  for  layer  2.  Layers  3  and  4  distributions  are 
quite  similar. 

Above  2000  m,  the  occurrence  of  severe  icing  is  greatly 
diminished.  Layer  6  (  250  1-3000  m)  shown  in  Figure  2.3,  is  repre¬ 
sentative  of  the  upper  four  layers.  The  occurrence  of  higher 
probabilities  is  quite  spotty  with  only  small  areas  in  the  north¬ 
ern  Soviet  Union  and  Norway  noteworthy. 

The  severe  icing  environment  is  of  very  limited  extent 
during  the  summer  as  shown  in  Figure  2.4  for  the  eighth  layer 
(450  1  -  6000  m).  This  is  the  only  summer  layer  in  which  high 
icing  probabilities  appear.  The  exceptions  are  for  the  South 
China  Sea  and  over  the  northern  Bay  of  Bengal  and  over  central 
Scandinavia.  The  high  probabilities  associated  with  the  Aslan 
monsoon  are  probably  spurious  as  the  3DNEPH  analyses  in  this 
region  overestimates  cloud  amount  (Hughes,  Henderson-Sel lers 
1985)  . 

The  severe  icing  environment  during  autumn  shown  in  Figures 
2.5,  2.6,  and  2.7  is  located  even  more  northward  than  its  spring 
counterpart.  The  most  extensive  occurrences  shown  in  Figure  2,5 
for  the  501  -  1000  m  layer  are  located  over  Alaska,  northern 
Quebec  and  Baffin  Island,  northern  Norway,  and  the  northern  Soviet 
Union.  The  distribution  for  the  fourth  layer  (1501  -  2000  m) 
shown  in  Figure  2.6  is  quite  similar,  howe  ■■■;.,  the  severe  icing 
environment  is  more  extensive  over  the  north  A  ’antic  from 
Greenland  to  Norway  and  over  eastern  Siberia  and  the  Kamchatka 
peninsula.  Layers  3  and  4  distributions  are  again  quite  similar. 
Figure  2.7  shows  the  distribution  in  layer  6  which  is  again  quite 
representative  of  the  layers  above  2000  m.  Higher  probabilities 
occur  over  Alaska  and  the  northern  Soviet  Union. 

The  northern  hemisphere  severe  icing  environment  is  most 
widespread  during  the  winter  as  shown  in  Figures  2.8,  2.9,  and 
2.10.  Figure  2.8  shows  the  distribution  for  layer  2  (50  1  -  1000 
m)  for  which  there  are  three  large  areas  of  severe  icing:  the 
eastern  Asia  littoral  from  Japan  to  the  Kamchatka  peninsula, 
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Occurrence  of  Supercooled  Liquid  Water  Content  Greater 
than  1gm/m*  in  April  at  2501-3000  meters.  (Map  shows 
only  those  regions  where  the  frequency  of  occurrence  is 
101  or  greater.) 
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F igure 


Occurrence  of  Supercooled  Liquid  Water  Content  Greater 
than  1gm/m*  in  July  at  4501-6000  meters.  (Map  shows 
only  those  regions  where  the  frequency  of  occurrence  is 
10$  or  greater  . ) 


Occurrence  of  Supercooled  Liquid  Water  Content  Greater 
than  1gm/m’  in  October  at  1501-2000  nr.eters.  (Map  shows 
only  those  regions  where  the  frequency  of  occurrence  is 
’•OS  or  greater.) 
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Figure  2.7  Occurrence  of  Supercooled  Liquid  Water  Content  Greater 
than  Igm/m’  in  October  at  2501-3000  meters.  (Map  shows 
only  those  regions  where  the  frequency  of  occurrence  is 
10%  or  greater  ,  ) 
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Figure  2.8  Occurrence  of  Supercooled  Liquid  Water  Content  Greater 
than  Igm/m1  in  January  at  501-1000  meters.  (Map  snows 
only  those  regions  where  the  frequency  of  occurrence  is 
10%  or  greater . ) 
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Figure  2.9  Occurrence  of  Supercooled  Liquid  Water  Content  Greater 
than  Igm/m*  in  January  at  1501-2000  meters.  (Map  shows 
only  those  regions  where  the  frequency  of  occurrence  is 
10%  or  greater . ) 
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Figure  2.10  Occurrence  of  Supercooled  Liquid  Water  Content  Greater 
than  Igm/m*  in  January  at  2501-3000  meters.  (Map  shows 
only  those  regions  where  the  frequency  of  occurrence  is 
10%  or  greater.) 


across  the  central  United  States  and  eastern  Canada  to  the  western 
north  Atlantic,  and  the  eastern  north  Atlantic  from  Greenland  to 
Scandinavia  and  the  North  Sea.  Smaller  areas  are  located  over  t ^ 
northern  Soviet  Union,  Alaska,  and  western  Canada.  The  winter 
severe  icing  environment  seems  to  coincide  with  the  northern  hemi¬ 
sphere  cyclone  tracks  and  is  highly  correlated  with  the  higher 
probabilities  of  cloud  as  revealed  the  Hughes  and  Henderson- 
Sellars  cloud  climatologies. 

It  might  be  surprising  that  the  severe  icing  environment  is 
of  such  limited  extent  over  most  of  the  Soviet  Union.  However, 
extensive  moisture  sources  are  located  far  from  the  Soviet  Union 
and  the  low  temperatures  there  additionally  limit  moisture.  This 
situation  is  apparently  repeated  for  much  of  western  and  central 
Canada  where  severe  icing  probabilities  are  low. 

The  distribution  for  layer  4  (1501  -  2000  m)  shown  in  Figure 
2.9  is  similar  to  that  of  Figure  2.8  and  similar  to  that  for  layer 
3  as  well.  Overall  probabilities  are  higher  and  the  severe  icing 
environment  over  central  and  southeastern  Europe  is  considerably 
more  extensive.  Some  of  these  latter  locations  have  altitudes  of 
1000  m  or  greater  so  that  layer  4  may  be  indicative  of  conditions 
not  far  above  cloud  base  for  this  region. 

Finally,  Figure  2.10  shows  the  winter  distribution  for  layer 
6  (  250  1  -  3  000  ra).  Again  this  layer  is  representative  of  the  up¬ 
per  4  layers.  Severe  icing  is  limited  to  eastern  Greenland,  the 
northeastern  Norway  coast  and  a  few  spots  in  Europe  and  Asia.  As 
before,  the  probability  of  severe  icing  occurrence  decreases  with 
increasing  altitude  so  that  the  probability  is  virtually  nil  for 
layer  eight. 

As  mentioned  before,  the  frequencies  of  severe  icing  occur¬ 
rences  calculated  using  the  ETAC  LWC  data  base  seem  high  and  this 
is  attributed  to  use  of  the  original  Sraith-Feddes  model. 
Nonetheless,  the  distribution  of  severe  icing  especially  for 
winter  coincides  with  high  frequency  of  cloud  occurrence  as 
revealed  in  the  Hughes  and  Henderson-Sellars  cloud  climatology. 

We  believe  that  in  spite  of  the  caveats  above,  the  severe  icing 
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climatology  as  described  here  has  validity  in  as  far  as  the  dis¬ 
tribution  of  severe  icing  is  concerned.  Thus,  we  can  assess  how 
representative  are  the  existing  FAA/NRL  icing  data  base  and  for¬ 
eign  data  additions  thereto  are  of  the  northern  hemisphere  severe 
icing  environment. 

2.4  AREAS  OF  THE  NORTHERN  HEMISPHERE  SEVERE  ICING  ENVIRONMENT 

WITH  LTTTLE  OR  NO  DATA  COVERAGE 

Figure  2.11  shows  the  locations  of  data  collection  for  the 
current  FAA/NRL  icing  data  base  and  the  foreign  data  colleoted  un¬ 
der  the  current  effort.  The  data  has  been  mainly  collected  in 
Europe  and  North  America  but  also  incluaes  data  collected  in  South 
Africa.  Extensive  data  collection  over  Lake  Miohlgan,  Nova  Scotia 
and  the  Syracuse,  New  York  area  seem  appropriate  as  the  frequency 
of  severe  icing  is  relatively  high  for  these  locations.  On  the 
other  hand,  the  highest  probabilities  of  severe  icing  are  looated 
over  Lake  Superior,  the  Saint  Lawrence  river  valley,  British 
Columbia,  the  North  Saa,  and  In  proximity  to  the  Japanese  and 
Sakhalin  islands.  Data  for  these  particular  locations  are  few  or 
nil.  Thus,  there  is  the  vexing  question  as  to  whether  an  extreme 
value  analysis  of  all  existing  data  might  inadequately  portray 
severe  icing.  It  is  obviously  impossible  to  make  measurements  at 
all  locations  where  the  maps  in  this  seotlon  suggest  a  high  prob¬ 
ability  of  LWC.  But  the  maps  do  suggest  that  additional  data 
collection  efforts  for  Lake  Superior,  the  Saint  Lawrenoe  river 
valley,  the  North  Atlantic  and  North  Sea,  and  the  Aslan  littoral 
would  probably  yield  Important  information.  Because  of  the  high 
coat  of  aircraft  measurements  the  alternate  data  development 
scheme  discussed  in  Seotlon  3  in  conjunction  with  RTNEPH  and 
northern  hemisphere  analysis  fields  should  be  implemented.  In  ad¬ 
dition,  satellite  borne  and  ground-based  microwave  sounders  could 
prove  useful  In  calibrating  the  Smlth-Feddes  model. 

It  is  also  possible  to  follow  Jeok's  (1983)  suggestion  of 
deciding  if  there  is  sufficient  data  by  examining  the  maximum  LWC 
measured  for  a  particular  v-eather  condition  and  then  letting  that 
guide  the  answer  to  the  question  "How  much  data  is  enough?" 
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This  i s  at  present  an  impossible  question  to  answer  for  the 
severe  icing  areas  with  no  measurements  since  there  is  no  maximum 
LWC  to  guide  the  answer.  The  modified  updated  Smlth-Feddes  model 
could  help  answer  such  a  question  by  predicting  the  highest  con¬ 
ceivable  LWC.  A  principle  difference  between  the  foreign  data  and 
the  existing  data  base  are  the  high  LWC's  measured  in  South 
African  convective  clouds  and  the  frequent  occurrence  of  LWC  >  1 

■3  0 

gra/nr  in  some  areas.  The  maximum  appears  to  be  3 .0  gra/mJ  with  two 
questionable  values  of  5.2  gm/m^.  Similar  values  have  been 
reported  by  the  Canadians.  According  to  the  Jeck  criterion  for 
data  miles: 


M  -  t>(LWCMAX) 


where 


2 

b  -  500  nmi/gm 

•j 

A  LWCmax  of  3.0  gm/nr  translates  to  1,500  nmi  of  data  col¬ 
lection,  Obviously  that  amount  of  data  collection  in  deep 
convective  clouds  is  expensive.  The  alternate  data  development 
discussed  in  Section  3  is  perhaps  more  realistic  (and  certainly 
less  expensive)  especially  for  deep  convective  clouds.  At  present 
it  is  impossible  to  assess  the  amount  of  data  collected  for  the 
various  weather  types  mentioned  by  Jeck  in  his  Table  B-  1  (  1983). 

It  seems  probable,  though,  that  all  the  weather  types  Jeck  indi¬ 
cates  as  deficient  in  measurement  will  continue  to  remain  so  even 
with  the  newly  collected  data.  The  difficulty  of  remedying  this 
deficiency  is  acute  for  these  weather  conditions  as  they  are  lo¬ 
calized  and  changeable  and  thus  hard  to  encounter.  Thus  the 
alternate  data  development  method  in  the  following  section  is  an 
attractive  means  for  establishing  not  only  how  much  data  is  needed 
but  also  to  extend  the  validity  of  limited  measurements  to  satisfy 
these  data  requirements. 
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SECTION  3 

ALTERNATE  DATA  DEVELOPMENT  SCHEME 


Current  aircraft  icing  research  is  concerned  with  ice 
protection  systems  that  are  energy  efficient  and  the  needs  of 
relaxing  the  certification  standards,  an  improved  characterization 
of  the  worldwide  icing  environment,  and  improved  forecasting  of 
icing  conditions  is  required.  The  worldwide  characterization  cur¬ 
rently  consists  of  retrieval  of  data  in  which  the  icing  parameters 
have  been  measured  in  addition  to  other  information  such  as  type 
of  cloud,  horizontal  or  vertical  extent  of  clouds,  and  synoptic 
weather  conditions.  As  discussed  in  Section  2,  the  parameter 
variation  and  especially  geographical  coverage  of  these  measure¬ 
ments  is  Insufficient  to  guarantee  that  the  resulting  data  base  is 
at  least  conservative.  To  guarantee  this  and  to  increase  the 
geographical  validity  either  additional  flight  data  surveys  or  an 
alternate  means  for  satisfying  the  data  requirements  are  neces¬ 
sary  . 

One  of  the  issues  of  the  original  contract  was  to  determined 
if  there  was  a  need  for  further  icing  data  measurements.  The  maps 
in  Section  2  show  that  there  are  extensive  geographical  areas 
where  severe  icing  is  expected  and  for  which  there  are  few  or  no 
measurements.  In  addition,  following  the  lead  of  Jeck  (1983)  cer¬ 
tain  types  of  meteorological  conditions  such  as  warm  fronts, 
occluded  fronts,  and  lake  effect  clouds  are  insufficiently  repre¬ 
sented  by  the  data  base.  Thus,  a  need  to  obtain  additional  data 
to  extend  the  data  base's  geographical  and  weather  condition 
validity  seems  justified.  The  FAA  anticipated  this  situation  by 
stipulating  that  if  additional  data  is  required,  an  alternate  data 
development  scheme  be  formulated  to  reduce  or  negate  the  reliance 
on  dedicated  and  costly  flight  surveys.  The  University  of  Dayton 
initiated  a  review  of  candidates  for  the  alternate  data  develop¬ 
ment  scheme.  The  results  of  this  review  are  included  in  Sections 
3 .  2  and  3.3* 

Based  on  this  review,  the  updated  Smith-Feddes  model  (SFM) 
was  selected  as  being  the  best  candidate  for  an  alternate  data 


development  scheme  baaed  on  a  validation  using  the  present  data 
base  of  LWC  and  MVD  values.  Section  3 .  **  presents  the  results  of 
the  validation.  Finally  Section  3.5  presents  two  strategies  by 
which  a  corrected  SFM  can  satisfy  the  requirements  of  an  alternate 
data  development  thus  negating  or  certainly  reducing  the  need  for 
a  number  of  costly  flight  surveys. 

3.1  THE  ICING  PARAMETERS 

Aircraft  icing  occurs  when  an  aircraft  encounters  super¬ 
cooled  liquid  water.  This  most  frequently  occurs  in  clouds  but 
can  also  occur  for  the  extremely  dangerous  situations  of  freezing 
drizzle  or  rain.  In  either  case  the  supercooled  liquids  freeze 
upon  impacting  the  aircraft  and  the  resulting  accretion  can 
seriously  degrade  the  aircraft  aerodynamio  performance,  impair  the 
operation  of  aircraft  instruments,  and  reduce  the  aircraft's  en¬ 
gine  performance.  The  amount  of  performance  degradation  is 
dependent  both  upon  the  type  of  aircraft  and  environmental  condi¬ 
tions.  For  a  given  aircraft  three  critical  environmental 
parameters  have  been  identified  (Pass,  198*0!  air  temperature 
(T),  liquid  water  content  (LWC),  and  drop  size  distribution  (DS). 
In  addition,  it  is  necessary  to  know  the  cloud  water  phase,  P^, 
because  relatively  large  areas  of  cloud  may  bo  glaciated  present¬ 
ing  a  benign  icing  environment. 

Temperature  is  Important  because  the  type  of  loe  accretion 
changes  dramatically  between  temperatures  near  freezing  and  those 
considerably  below  freezing.  At  temperatures  near  freezing  the 
drops  do  not  freeze  instantaneously  but  run  back  slightly  forming 
a  dear  glaze  ice  that  can  attain  spectacular  formations.  For 
temperatures  well  below  freezing  the  drops  freeze  instantaneously 
forming  a  rime  ice  which  is  not  quite  as  aerodynamlcally 
troublesome  as  the  glaze  ice.  Additionally,  for  colder  tempera¬ 
tures  the  available  LWC  is  normally  less  than  for  near  freezing 
temperatures. 

Of  the  threo  parameters,  DS  presents  the  greatest  difficulty 
in  both  measurement  and  calculation.  Drop  size  distributions  are 
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sensitive  to  the  cloud  condensation  nuclei,  strength  of  vertical 
motion,  mixing  of  clouds  and  cloud  free  air,  and  cloud  lifetime. 
These  factors  have  defied  the  concerted  attention  of  a  number  of 
investigators  so  that  there  is  no  general  method  for  their  ac¬ 
counting.  In  addition,  it  is  currently  difficult  or  impossible  to 
specify  vertical  motion,  cloud  mixing,  or  cloud  lifetime. 

Fortunately  there  may  be  a  way  out  of  this  dilemma.  Newton 
(1978  )  presented  results  that  support  the  substitution  of  a  median 
volume  drop  (MVD)  in  place  of  a  drop  spectrum.  Theoretical  and 
experimental  studies  can  determine  the  reasonableness  of  such  an 
approximation.  With  an  assumed  drop  size  it  is  still  necessary  to 
measure  or  infer  the  temperature  and  the  cloud  phase.  LWC  is 
critically  Important  as  it  is  directly  related  to  the  degree  of 
ice  accretion,  and  phase  determines  whether  the  cloud  particles 
will  adhere  to  the  aircraft. 

3.2  MODELING  CANDIDATES 

3.2.1  Smlth-Feddes  Model 

The  Smith-Feddes  Model  predicts  estimates  of  LWC, 
cloud  phase,  and  drop  size  distribution.  The  output  grid  is  the 
same  as  the  Air  Force  Global  Weather  Central's  (AFGWC)  3DNEPH 
(three-dimensional  nephanalysis )  with  15  vertical  layers  whose 
thickness  varies  with  altitude.  The  model  input  consists  of  the 
3  DNEPH  analysis  of  cloud  type,  the  AFGWC  analysis  of  grid  tempera¬ 
tures  interpolated  to  the  3DNEPH  grid,  as  well  as  information 
about  preolpltation.  For  a  given  cloud  type  and  temperature  there 
is  a  maximum  LWC  which  is  tabularly  assigned  by  the  model.  This 
maximum  LWC  is  modified  by  the  in  cloud  location  of  the  grid  point 
and  the  oloud  type. 

The  portion  of  cloud  water  which  is  unfrozen  varies 
from  1 00%  at  0°C  to  0%  at  -4Q°C  and  is  based  on  Russian  data. 
Finally,  drop  size  distributions  are  assigned  according  to  the 
oloud  type  following  the  distributions  presented  by  Diem  (19*18). 
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3.2.2  Updated  Smlth-Feddes  Model 

Rogers  et  al.  (1985a)  detail  the  updating  of  the 
original  Sraith-Feddes  model  discussed  above.  The  updated  version 
accepts  input  from  the  improved  cloud  analysis  model  ( RTNEPH )  and 
considerably  differs  in  microphysical  parameterizations  from  the 
original  model.  This  updated  version  is  based  on  the  most  recent 
microphysical  observations  available  at  the  time  of  the  above 
report.  There  is  no  need  in  this  report  to  detail  the  conversion 
from  accepting  3DNEPH  to  accepting  RTNEPH  input,  however,  the 
microphysical  changes  are  briefly  described  below. 

In  examining  data  on  the  vertical  distribution  of 
LWC,  the  authors  above  noted  that  the  profiles  in  some  cases  were 
close  approximations  of  the  adiabatic  LWC  profile  and  if  not  were 
some  fraction  of  the  adiabatic  LWC.  In  general,  they  believed 
that  stratus  and  stratocumulus  clouds  closely  profile  the 
adiabatic  profile  while  cumulus  clouds  diverted  with  height  from 
the  adiabatic  profile.  The  authors  believed  that  a  function  of 
the  adiabatic  profile  with  height  would  account  for  the  dirainish- 
raent  of  the  ratio  of  actual  LWC  to  adiabatic  LWC. 

The  implemented  adiabatic  computation  of  LWC  con¬ 
sists  of  several  steps.  From  cloud  base  the  moist  adiabatic  lapse 
rate  is  used  to  calculate  the  temperature  in  100  m  increments 
above  cloud  base.  These  temperatures  define  the  satuaration  vapor 
pressures  at  100  ra  increments  through  the  cloud.  Finer  resolution 
in  the  calculation  of  saturation  vapor  pressure  is  possible,  but 
is  not  significant  in  affecting  the  cloud  LWC.  The  difference  in 
saturation  vapor  pressure  over  a  100  m  interval  defines  the  addi¬ 
tional  condensed  moisture,  that  is  accumulated  beginning  at  cloud 
base  and  continuing  to  the  cloud  top.  Near  cloud  top,  LWC  is 
diminished  somewhat  to  account  for  enhanced  mixing  of  unsaturated 
air  above  the  cloud. 

Reduction  of  adiabatic  LWC  due  to  entrainment  fol¬ 
lows  Warner  (1970).  The  adjustment  consists  of: 

LWC /LWC  -  az  ♦  b  (3.1) 

3 
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where 


LWC  -  adiabatic  LWC , 

u 

z  -  height  above  cloud  base,  and 
a,b  -  parameters  which  are  In  turn  assigned  according  to  z. 

This  reduction  is  applied  in  the  event  of 
stratocumulus ,  altocumulus,  cumulus,  and  cumulonimbus  clouds.  No 
reduction  for  entrainment  was  originally  applied  for  stratus  or 
altostratus . 

In  the  Smith-Feddes  model  the  drop  size  distribution 
is  computed  from  a  parameterized  equation  so  that  the  number  of 
drops  in  a  given  drop  diameter  interval  is  the  fraction  of  total 
LWC  in  that  drop  size  interval  divided  by  the  mass  of  a  single 
droplet  of  that  size.  The  fractional  LWC  follows  the  Diem  dis¬ 
tribution  curves  and  is  based  on  cloud  type.  The  parametric 
equation  is  designed  primarily  to  give  the  correct  mode  of  the 
drop  size  distribution  as  well  as  the  percent  of  LW:  at  the  mode 
as  a  function  of  cloud  type.  Unfortunately,  as  shown  later,  it 
appears  that  this  approach  is  overspecified  and  as  a  result,  the 
method  greatly  overestimates  the  number  of  larger  drops. 

3.2.3  Adiabatic  LWC 

Adiabatic  LWC  can  be  calculated  from  rawtnsonde, 
satellite,  or  profiler  soundings  of  temperature  and  humidity.  The 
method  involves  calculating  the  cloud  base  temperature  r.n d  pres' 
sure.  With  temperature  and  pressure  known  at  this  levil,  it  is 
possible  to  infer  LWC  for  higher  heights  based  on  adiabatic  parcel 
ascent  from  this  level.  This  method  does  not  address  f.ntrainmenV. 
of  dryer  environmental  air  nor  does  it  allow  for  precipitation. 
Observations  of  cloud  LWC  values  typically  are  less  than  adiabatic 
and  the  ratio  of  actual  LWC  to  adiabatic  LWC  decreases  t  rom  oiouc. 
base  to  cloud  top.  The  validity  of  this  method  is  also  limited  to 
convective  clouds. 
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An  extension  of  this  method  has  been  developed  by 
the  Air  Weather  Service  (1969).  The  AWS  method  involves  using  a 
cloud  model  with  a  fixed  rate  of  entrainment  to  calculate  the  LWC 
by  height.  The  method  pertains  to  stratiform  and  convective 
clouds  depending  upon  the  convective  stability.  For  stratiform 
cloud  which  must  be  inferred  from  temperature  and  dew-point 
spread,  the  calculated  LWC  is  halved.  Finally  a  median  volume 
drop  diameter  (MVD)  size  of  1  -U  ym  for  stratiform  clouds  and  17  ym 
for  convective  clouds  is  specified.  These  values  are  closely  con¬ 
sistent  with  those  reported  by  Jeck  (1983)  who  summarized  the 
available  continental  United  States  observations  of  temperature, 
LWC,  and  DS.  Jeck  shows  that  MVD  values  vary  from  about  5  um  to 
about  30  ym  for  convective  clouds.  A  desirable  improvement  to 
this  method  would  be  to  specify  the  MVD  size  according  to  environ¬ 
mental  conditions.  Newton  (1978)  reported  good  comparison  between 
estimated  LWC  and  ice  accretion  using  this  method  and  that  ac¬ 
tually  measured. 

3.  2.  4  afcrl  Flr3t  Generation  Model  ( AFCRL-I ) 

Cunningham  and  Pierce  (197*0  employs  a  decision  tree 
approach  in  calculating  LWC  and  precipitation  phase.  Temperatures 
of  -’,00C,  -I5°C,  and  0CC  designate  respectively  transitions  be¬ 
tween  ice  crystals,  small  snow,  large  snow,  and  rain.  The  LWC 
values  are  typical  values. 

Tha  model  input  Includes  a  sounding  in  addition  to 
the  current  weather  and  intensity  of  precipitation;  the  type  of 
low  cloud  gives  an  indication  of  presence  or  absence  of  convective 
activity. 

3*  2.  5  AP.'rtL  Second  Generation  Model  (AFCRL-II) 

The  second  generation  model  depends  on  preparation 
of  vertical  time  height  aoross  sections  to  which  LWC  values  are 
added  '.Cunningham  and  Fierce  197*1).  The  cross  sections  include: 

(1)  Surface  observation  of  weather,  precipitation  type  and 
amount,  sky  cover,  cloud  type  and  amount,  cloud  base  heights, 
winds,  temperatures,  dew-point,  and  pressure. 
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(2)  Temperature,  dew-point  temperature  depression,  and 
winds  for  significant  and  special  levels  in  the  vertical. 

Satellite  visual  and  infrared  photographs,  surface  synoptic  maps, 
and  3 DNEPH  analyses  are  also  used.  The  cross  sections  are 
analyzed  for  stable  and  conditionally  unstable  layers  to  which 
cloud  patterns  are  added.  Values  of  LWC  and  hydroraeteor  con¬ 
centrations  are  assigned  according  to  published  values  in  the 
literature.  In  this  form  the  model  is  unsuitable  for  computer 
programming . 

3*2.6  Meteorology  Research  Inc.  Models  (MRI) 

The  MRI  models  MRI-RH  and  MRI-VV  are  respectively 
for  relative  humidity  and  vertical  velocity  (Heyrasfield  and  Howard 
1975).  These  models  represent  an  attempt  to  develop  nomograms  for 
assigning  LWC  by  cloud  and  weather  type.  The  nomograms  are  formu¬ 
lated  on  the  basis  of  vertical  moisture  flux  between  atmospheric 
levels.  For  the  vertical  velocity  nomogram  it  is  necessary  to 
calculate  a  profile  of  vertical  velocioty. 

Pierce  et  al.  (  1975  )  compared  the  performance  of 
several  of  the  above  model  (not  including  the  updated  Smith-Feddes 
model)  that  prescribe  LWC  with  a  limited  number  of  observations. 
These  are  shown  in  Figures  3*  la  and  3*15.  The  comparisons  are 
not  terribly  encouraging,  however,  as  shown  later  the  updated 
Smith-Feddes  model  shows  skill  in  predicting  LWC, 

3  •  2.  7  Matveev 

Matveev  (1984)  presented  a  discussion  on  the  cal¬ 
culation  of 

ZU 

Q,  -  /  p5dz  (3*  2) 

6  Z 

L 

where  Q{  is  the  integrated  cloud  water  content  and  Z^  and  Zy  are 
the  cloud  base  and  cloud  top  heights,  respectively.  The  method 
involves,  using  a  simplified  version  of  the  equation  for  conserva¬ 
tion  of  moisture  to  which  equations  for  heat  content,  vertical 
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Figure  3.1a. 
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velocity,  and  the  turbulence  coefficient  are  added.  For  this 
method  the  fall  velocity  of  cloud  elements  was  assumed  independent 
of  height.  The  resulting  equation  set  was  solved  numerically  by 
the  Runge-Kutta  method.  A  temperature  profile  from  a  '"adiosonde 
or  the  like  may  be  used  in  place  of  the  heat  content  equation. 
Matveev  claims  good  comparison  of  calculated  by  this  method  to 
values  measured  by  aircraft.  has  a  strong  dependence  on  ground 
level  temperature  which  has  been  verified.  However,  other 
theoretical  relationships  with  the  vertical  velocity,  turbulent 
mixing,  and  drop  fall  velocities  and  turbulence  coefficients  which 
require  diagnosis  or  calculation  by  other  sophisticated  models 
such  as  those  discussed  next.  Finally,  is  an  integrated  rather 
than  a  point  value. 

3.2.8  Other  Models 

Recent  developments  in  cloud  and  mesoscale  modeling 
(Kierap  and  Wilhelm3on  1978,  Jun  and  Hsu  1984)  enable  grid  point 
calculation  of  LWC.  The  calculations  range  from  microphysical 
pararaeteri zations  such  as  Kessler  (1969)  or  Sundqvist  (1981)  to 
detailed  microphysical  calculations  (Silverman  and  Glass  1985). 

In  the  latter  case  the  resulting  drop  spectrum  as  well  a3  LWC  is 
calculated.  For  convective  clouds  the  model  present  a  feasible 
means  for  alternate  data  de . elopme nt  .  however,  their  computational 
and  input  data  requirements  are  quite  stringent.  In  this  sense, 
the  Silverman  ana  Glass  1- dimensional  model  is  attractive  as  it 
only  requires  a  single  input  sounding;  on  the  other  hand,  most  1- 
dimensional  models  overpredic*'  LWC  for  cases  of  environmental  wind 
shear  (Cotton  and  Tripoli  1978)  so  such  a  model  has  limited 
validity.  For  stratiform  clouds  it  is  necessary  to  have  detailed 
3-diraen3 tonal  calculation  so  that  vertical  motions  can  be  faith¬ 
fully  portrayed. 

3.3  REMOTE  SENSING 

Passive  and  active  microwave  sensors  are  promising  means  for 
obtaining  LWC  in  clouds  (Devault  and  Katsaros  1983,  Rauber  et  al . 
1982),  however,  they  are  currently  restricted  to  integrated  values 
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(Shepard  1985).  Current  satellite  instruments  are  limited  in  the 
number  of  available  microwave  channels  and  this  is  likely  to 
remain  the-  case  into  the  next  decade.  Thereafter,  with  more  chan¬ 
nels  available  it  may  be  possible  by  inversion  techniques  to 
derive  a  vertical  profile  of  LWC. 

Ground  based  methods  for  calculating  LWC  seem  to  be  further 
along  in  development.  Rauber  et  al.  (1982)  report  on  favorable 
comparisons  between  radiometer  measured  values  of  LWC  to  those 
measured  by  a  research  aircraft.  One  goal  of  a  field  study 
recently  completed  in  Louisiana  was  to  test  a  tomographic  tech¬ 
nique  which  can  yield  a  profile  of  LWC  (Knight  1985).  It  will  be 
some  time  before  the  data  from  this  study  can  be  fully  evaluated. 

Although  remote  sensing  methods  seem  to  offer  good  future 
prospects  for  regular  observation  of  LWC,  it  is  currently  possible 
to  only  calculate  the  vertically  integrated  amount.  This  amount 
can  serve  to  calibrate  a  method  such  as  the  upoated  Smith-Feddes 
which  provides  vertical  detail  of  the  LWC  distribution.  In  addi¬ 
tion  for  ground  based  systems  there  13  the  problem  of  deployment 
especially  over  oceanic  and  mountainous  regions  which  unfor¬ 
tunately  seem  to  be  prime  locations  for  aircraft  icing  occurrence. 

3.4  VALIDATION  OF  THE  UPDATED  SMITH-FEDDES  MODEL 

Based  on  the  review  of  modeling  candidates  in  Section  3.2 
and  of  remote  sensing  techniques  in  Section  3.3.  it  became  ap¬ 
parent  that  the  updated  Smith-Feddes  model  offered  the  best 
promise  as  an  alternate  data  development  scheme.  The  model  was 
untested,  however,  so  a  thorough  validation  was  in  order  before  an 
unfettered  recommendation  could  be  made.  The  FAA/NRL  icing  data 
base  was  an  obvious  means  for  the  validation  provided  that  input 
data  was  available  to  satisfy  the  Smith-Feddes  requirements.  For 
a  number  of  data  base  measurements,  supplementary  data  from  the 
measuring  aircraft  or  a  conjunctive  rawinsonde  sounding  was  avail¬ 
able.  After  some  manipulation  it  was  possible  to  calculate 
temperature,  LWC,  and  drop  size  distribution  with  the  Smith-Feddes 
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model  and  then  compare  the  predicted  value  to  the  measured  value 
in  the  data  base. 


3.4/1  Preparation 


Smi th-Feddes 


Data 


One  of  the  options  for  running  the  updated  Smith- 
Feddes  model  is  via  a  card  image  file  (Rogers  1985b).  The 
required  file  consists  of  the  followings  latitude  and  longitude 
for  the  calculation  position;  present  weather,  visibility  and  time 
for  the  calculation  point;  total  percent  cloud  cover  and  partial 
percent  cloud  cover  for  each  cloud  deck  (up  to  4  decks);  coded  in¬ 
put  for  the  type  of  cloud;  cloud  base  and  cloud  top  heights  for 
each  cloud  deck;  surface  elevation  for  the  calculation  position; 
heights  in  meters  of  the  standard  pressure  levels  of  1000,  850, 
700,  500,  400,  300,  200,  150,  and  100  mb;  and  temperatures  in  de¬ 
grees  Kelvin  at  the  standard  pressure  levels  and  at  the  surface. 

In  many  of  the  data  base  entries  only  the  measured 
or  estimated  cloud  base  heights,  type  of  precipitation  if  any, 
surface  elevation,  and  cloud  type  were  available.  Thus,  sup¬ 
plementary  sources  were  necessary  to  obtain  standard  and  special 
pressure  level  temperatures  and  heights.  In  some  cases,  the 
aircraft  measured  temperature,  pressure,  and  altitude  in  conjunc¬ 
tion  with  daily  weather  map  surface  and  500  mb  analyses  were 
sufficient.  In  the  case  of  the  Spanish  PEP  data  it  was  necessary 
to  obtain  the  standard  and  special  pressure  level  heights  by 
hydrostatically  integrating  the  available  temperature  and  pressure 
soundings,  ilegardless  of  these  difficulties  a  total  of  about  100 
comparisons  were  made.  This  represents  about  50$  of  the  total 
measurements  included  in  the  data  base.  It  was  not  possible  to 
use  any  NACA  measurements  since  supplementary  data  for  standard 
level  heights  and  temperatures  was  unavailable. 

3.4.2  Philosophy  for  Comparison  of  Smlth-Feddes  Calculated 

LWC,  T,  and  PS  to  FAA/NRL  Observations 

The  FAA/NRL  data  base  of  supercooled  LWC  measure¬ 
ments  is  organized  according  to  the  collecting  agency,  time, 
cloud,  and  event.  A  given  "cloud”  corresponds  to  one  cloud  base 
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height  and  environmental  sounding  and  might  include  several 
"events".  (See  Jeck  1983  fcr  details  of  cloud  and  event  defini¬ 
tion.)  Since  the  Smith-Fedd^s  model  was  originally  designed  for 
data  with  25  nautical  mile  horizontal  resolution,  we  decided  to 
compare  average  FAA/NRL  "cloud"  values  of  LWC,  T,  and  DS  to  the 
Smith-Feddes  calculations.  For  a  cloud,  there  may.be  several 
events  which  are  representative  of  a  finer  aoale  than  is  ap¬ 
propriate  to  the  Smith-Feddes  calculations.  In  addition,  this 
larger  scale  is  more  relevant  to  aircraft  icing  since  an  airoraft 
is  not  immediately  disabled  by  penetrating  small  regions  of  high 
LWC  content  nor  is  it  immediately  unburdened  of  loe  when  transit¬ 
ing  small  scale  unsaturated  regions.  The  general  trond  of  loe 
accretion  over  a  10  mile  or  larger  distance  should  then  oorrelate 
reasonably  well  with  average  cloud  values  as  described  above. 

3.^.3  Changes  to  Updated  Smith-Feddes  Modol 

During  an  Initial  validation,  several  deficiencies 
of  the  updated  Smith-Feddes  wore  noted.  The  model  appeared  to 
calculate  slightly  too  muoh  LWC  for  stratus  olouds  and  too  little 
LWC  for  stratocumulus  and  cumulus  clouds.  In  addition,  con¬ 
siderable  scatter  between  calculated  and  observed  LWC  values  was 
noted  for  all  clouds.  Smith-Feddes  calculated  MVD  values  were  in 
significant  disagreement  with  those  observed  and  the  Smith-Feddes 
drop  size  calculation  procedure  was  identified  as  unlikely  to 
produce  realiatlo  MVD  values. 

The  LWC  discrepancies  were  attributed  to  two 
sources:  inoorreot  cloud  base  temperature  or  improper  aooountlng 

for  the  mixing  of  cloudy  and  unsaturated  air.  The  first  source 
was  ameliorated  by  allowing  the  model  greater  vortical  resolution 
of  input  temperature  and  heights  versus  pressure.  This  allowed  a 
properly  Identified  oloud  base  temperatures  whloh  are  associated 
with  frontal  inversions  or  the  like  which  are  poorly  resolved  by 
standard  pressure  lovel  data  alone. 

The  second  source  was  not  as  easily  handled,  In 
general,  sub-adlabatlo  LWC  values  exist  in  the  bulk  of  oonveotive 
clouds  and  also  apparently  for  many  layer  clouds.  Precipitation 
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and  entrainment  of  dryor  external  air  oontrlbutea  to  the  diralnish- 
ment  of  LWC  values.  Entrainment  may  take  the  form  of  organized 
Inflow  of  external  air  similar  to  that  observed  for  laboratory 
plumes,  via  turbulent  mixing  driven  by  shear  and  buoyant  energy 
production,  or  via  internal  penetrative  plumes  whloh  arise  from 
mixing  near  the  oloud  top  of  cloudy  and  oloud  free  air.  Preoise 
modeling  of  ail  these  prooesses  is  currently  impossible  even  in 
the  very  sophi eticcted  three-dimensional  models  of  deep  oonveotion 
(see  Section  1.2. Hi.  Crude  but  reasonable  parameterization  of 
oloud  mlcrophysios  and  turbulenoe  is  necessary.  Favorable  com¬ 
parisons  between  observed  and  calculated  LWC  values  for  layer 
clouds  were  obtained  for  the  Smith  Feddes  model  when  the  Warner 
reduction  curve  was  used  for  stratus  olouds  and  the  Skatskil 
reduotlon  curve  for  atratooumulua  clouds  (Figure  3.2).  For  oon- 
veotlve  clouds  the  Skatskil  (1965)  reduotlon  ourve  gives  superior 
results  to  the  Warner  curve.  Despite  the  reasonable  LWC  valuos, 
the  calculated  drop  size  distributions  remain  uniformly  unrealis¬ 
tic. 

We  have  not  currently  rectified  the  drop  distribu¬ 
tion  problem,  however,  we  can  suggest  a  procedure  that  is  expected 
to  greatly  Improve  Da  calculation.  Following  Johnson  (1986)  and 
Berry  and  Kelnhart  (1974)  the  oloud  droplet  distribution  in  many 
oases  can  be  described  by  a  gamma  distribution  as 

f<x>  •  ^  GOO  aY  e"(1*Y)3  (3.3) 


where 


S 


Nx 

UW 


GH) 


(3.4) 

(3.5) 


and  where  f<x>  Is  the  number  density  of  drops  of  mass  x,  N  is  the 
total  number  of  oloud  droplets  per  unit  volume  and  7  is  the  shape 
parameter  which  specifies  the  distribution  breedth.  Although  7-2 
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Figure  3.2  Ratio  of  Observed  Liquid  Water  Content  (A) 
to  Adiabatic  Liquid  Water  Content  (Q  )  as  A 
Function  of  Height  Above  Cloud  Base  from 
Warren  (1970) . 


seems  appropriate  for  many  observed  distributions,  it  might  be 
necessary  to  reduce  Y  with  height  down  to  as  low  as  0  to  account 
for  the  behavior  of  some  observed  distributions  with  height.  It 
may  be  possible  to  relate  Y  to  environmental  parameters  such  as 
cloud  location,  height  of  cloud  base,  season,  wind  shear,  CCN , 
etc . 

At  any  rate,  integrating  equation  (3.3)  over  x  can 
yield  the  total  cloud  droplet  number  per  unit  volume 

N  -  /  f < x>  dx  (3.6) 

and  the  liquid  water  content 

LWC  -  /  x  f<x>  dx.  (3.7) 

With  LWC  already  known  from  the  Sraith-Feddes  calculations  it  is 
possible  to  infer  N  and  hence  the  drop  size  MVD.  The  key  issue  is 
proper  selection  of  Y  in  terras  of  environmental  parameters, 
however,  even  a  default  calculation  in  which  Y  is  specified  as  2 
would  yield  useful  results. 

3.^.4  Comparison  of  Smlth-Feddes  Calculated  LWC,  T,  and  PS 

Values  to  rAA/hlftL  Obser'vatio~ns 

Table  3.1  lists  the  calculated  and  observed  tempera¬ 
tures,  MVD  and  LWC  values  for  all  validation  cases.  The  Table  3.1 
arrangement  corresponds  to  the  order  in  the  FAA/NRL  data  base  for 
ease  in  cross  referencing.  Each  entry  in  Table  3.1  lists  the 
location  and  collecting  agency  of  the  observation,  the  cloud  num¬ 
ber  which  corresponds  to  the  FAA/NRL  data  base  entry,  the  date, 
cloud  type,  average  LWC  measured  respectively  by  the  Johnson- 
Williaras  and  FSSP  devices,  the  LWC  calculated  by  the  model,  the 
mean  altitude  of  the  observation,  and  the  calculated  and  observed 
cloud  base  temperature.  In  general,  these  latter  two  numbers 
agree  closely,  however,  in  a  few  instances  it  was  difficult  for 
the  model  to  properly  resolve  strong  inversions.  Cloud  base  tem¬ 
perature  difference,  however,  i3  not  the  main  contributor  to 
differences  in  observed  and  calculated  LWC  values.  As  mentioned 


3-15 


COMPARISON  OF  CALCULATED  AND  OBSERVED  TEMPERATURE.  MVD 
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TABLE  3.1  (continued) 
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TABLE  3.1  (concluded) 


©• 


s  * 

m 

*1 

0 

S' 

O 

m 

o 

0 

O 

# 

> 


3-26 


previously,  the  strength  and  nature  of  cloud  nixing  varies  con¬ 
siderably  depending  on  the  environment  and  the  Warner  and  Skatskii 
reduction  curves  are  only  mean  summaries  of  the  reduction  which 
has  an  obvious  varlanoe.  In  addition,  there  is  instrument  error 
whloh  is  evident  when  the  two  looated  measuring  devices  are  com¬ 
pared.  In  general,  the  LWC  comparisons  in  Table  3.1  for  a  variety 
of  environments,  collecting  agencies,  and  cloud  types  are  quite 
good . 

In  our  investigation  we  looked  at  several  methods  of 
evaluating  the  comparisons  of  the  Smlth-Feddes  predicted  LWC 
values  to  those  observed  by  the  J ohnson* Wll 1  lams  and  partlole 
spectrometer  measuring  devices.  We  initially  performed  a  least 
squares  fit,  on  the  raw  data,  to  the  line  of  the  form 

Y  •  m  X  •  b .  (3.8) 

Th*  rmmulim  from  thl«  method  gay*  y§  an  indication  of  any  Dias 
that  the  model  might  have.  For  the  oases  whloh  compared  the 
Smlth-Feddes  predicted  LWC  to  the  Johnson-Willlams  and  partlole 
spectrometer  measured  LWC  (Figures  3*3  and  3.4,  respectively )  for 
all  cloud  types  considered  (St,  3c,  Ao,  and  Cu),  the  linear 
regressions  produced  the  equations, 

3'F  •  0.05  (J-W)  *  0.08  (3.9) 

and 

3-r  -0,01  (P-S)  •  0.09  .  (3.  10) 

Thus,  we  see  that  the  blse,  If  any,  In  ths  model  is  lrislgnlf  losnt . 
For  this  reason,  snd  the  feet  thst  if  the  measuring  devices  did 
not  register  any  LWC  at  all  thers  would  bo  no  clouds  and  If  there 
are  no  clouds,  then  the  Smlth-Feddes  model  would  also  Indicate  no 
LWC,  Thus,  we  force  our  regression  line  through  the  origin.  The 
form  of  the  regression  line  then  beoomes, 

Y  •  m  X  .  (3.11) 
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Measured  LWC  vs.  S~F  Predicted  LWC  for  ST,  SC,  AC 
CU  Cloud  Types. 


Heasured  LWC  vs 
CU  Cloud  Types. 


Figures  3.3  and  3.4  show  that  for  larger  values  of 
liquid  water  oontent,  the  data  points  become  more  spread  cut. 

This  spreading  out,  or  fanning,  of  the  data  indicates  that,  in 
this  form,  the  data  does  not  lend  Itself  to  a  regression  analysis, 
because  the  more  extreme  points  have  a  larger  Influence  on  the 
regression  line.  However,  if  we  go  to  a  logarithmio  soale,  as  in 
Figure  3.5,  we  see  that  the  spreading  of  the  data  is  more  uniform, 
and  therefore  the  influence  of  the  points  will  be  more  even. 

Thus,  the  form  of  the  regression  line  transforms 
from  Eq.  (3.11)  to 

tn  X  -  in  X  ♦  tn<m )  ,  (  3.  12) 

but  since  m  is  a  oonatant,  tn(m)  is  a  constant,  so  let  ln(m)  -  b. 
Eq.  (3>12)  can  then  be  rewritten  as, 

tn(3-F )  -  tn(O-W)  ♦  b.  (3.  13) 

From  here  we  oan  solve  directly  for  the  intoroept,  b, 

b  •  .  (  3.  14) 

The  linear  regression  oaloulation  for  eaoh  of  the  comparison  oases 
reduces  to  finding  the  mean  for  b,  (b).  Thus, 

in(S-F)  -  m(J-W)  *5  (3.  15) 


is  the  regression  line  in  the  in  vs.  in  space.  To  return  to  the 
unsealed  space,  all  we  need  to  do  is  perform  the  inverse  logarithm 
on  both  sides  of  Eq.  (3.14)  to  obtain, 


where 


3 ~V  -  #b  J -W  or  li-r  •  a  J-W  , 


m 


(3.  10) 


(3.  17  ) 
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This  method  was  used  to  calculate  the  regression  lines  for  each  of 
the  comparison  cases  and  the  results  are  displayed  in  each  of 
their  corresponding  Figures  (.  ..  -  3.14). 

After  obtaining  the  regression  lines,  the  present 
error  in  the  Smith-Feddes  model  could  be  determined  for  each  of 
the  comparison  cases.  The  percent  error  for  each  case  is  the  mean 
of  the  percent  error  of  each  of  its  individual  data  points.  Thus, 
if  Y  is  the  actual  Smith-Feddes  value  and  Y'  is  its  corresponding 
predicted  value  (Y'  -  m  X,  where  X  is  the  corresponding  Johnson- 
Williams  or  particle  spectrometer  measured  LWC),  and  n  is  the 
number  of  data  points  in  the  comparison  case,  then  the  percent  er¬ 
ror  is  calculated  by, 

inn  ^  I  y  *  -  y  i 

Percent  Error  -  — — -  £  -* — - — L  .  (3.  18) 

n  i- 1  Y 

The  percent  error  calculations  for  each  of  the  comparison  cases 
gives  an  indication  of  how  well  the  S-F  model  predicts  the  LWC. 

Figure  3.3  shows  the  particle  spectrometer  measured 
LWC  versus  the  Smith-Feddes  calculated  LWC  for  all  clouds.  Most 
of  the  values  1J*3  i<thin  the  dashed  line.  However,  those  values 
that  do  not,  ind.^ce  the  Smith-Feddes  overpredicts  LWC.  Figure 
3.4  shows  the  same  comparison  except  for  the  Johnson-Will iams 
device.  The  results  in  Figures  3-3  and  3.4  are  quite  similar  with 
the  Smith-Feddes  agreeing  slightly  better  with  the  particle 
spectrometer  (PS)  than  with  the  Johnson-Williams  (JW).  The  per¬ 
cent  errors  between  model  and  observation  are  nearly  the  same  as 
for  the  lnt6rcoraparison  of  measuring  devices. 

Figures  3.6  and  3.7  show  the  LWC  comparisons  for 
respectively  the  PS  and  the  JW  device  for  stratus  clouds.  The 
percent  errors  are  slightly  lower  than  for  all  the  clouds  com¬ 
parison  however,  the  largest  stratus  LWC  is  only  about  0.2gm/m^  so 
that  these  percent  errors  translate  to  fairly  low  absolute  errors. 
Some  measurement  errors  may  be  independent  of  LWC  and  are  rela¬ 
tively  large  for  low  LWC  values.  As  is  the  case  for  the  all  cloud 
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Measured  LWC  vs.  S-F  Predicted  LWC  for  Stratus 


comparison,  the  Smith-Feddes  calculated  values  agree  slightly  bet¬ 
ter  with  the  particle  spectrometer. 

Figures  3.8  and  3.9  show  the  LWC  comparisons  for 
the  PS  and  the  JW  for  stratocumulus  (Sc)  clouds.  The  model  con¬ 
tinues  to  agree  more  closely  to  the  particle  spectrometer  and  for 
Sc  the  agreement  is  particularly  good  even  for  higher  LWC  values. 

Most  of  the  comparison  values  are  for  observed  values  below 

3 

0.3gm/ra  so  it  would  be  useful  to  have  additional  Sc  measurements 

3 

where  LWC  was  greater  than  0.3gm/m  . 

Figures  3*10  and  3-11  show  the  LWC  comparisons  for 
the  PS  and  the  JW  for  cumulus  (Cu)  clouds.  The  percent  errors  are 
intermediate  between  those  for  stratus  and  those  for  stratocumulus 
and  the  comparison  between  the  JW  and  the  model  calculation  is 
slightly  better  than  for  the  PS.  Most  of  the  LWC  comparison  which 
lie  significantly  outside  the  percent  error  lines  occur  for 
cumulus.  There  is  an  additional  factor  for  cumulus  which  has  not 
been  previously  mentioned.  In  environmental  wind  shear  new  growth 
of  cumulus  clouds  occurs  on  the  upshear  side  of  the  cloud,  while 
the  downshear  side  consists  mainly  of  dissipating  cloud  (Rogers  et 
al.  1985).  As  a  result  the  maximum  LWC  values  occur  on  the  upwind 
side  of  the  cloud  and  gradually  diminish  downwind.  For  more  com¬ 
plex  shear  scenarios  the  maximum  LWC  center  will  be  shifted 
offcenter.  An  aircraft  may  not  adequately  sample  a  cumulus  formed 
in  wind  shear  because  of  the  direction  in  which  the  aircraft 
penetrates  the  cloud.  An  along  shear  path  is  more  likely  to 
faithfully  portray  the  cloud  LWC  as  opposed  to  an  across  shear 
path  which  has  little  chance  of  giving  the  true  cloud  averaged 
LWC.  Cumullform  LWC  structures  and  aircraft  penetration  trajec¬ 
tories  would  appear  to  be  a  factor  influencing  the  model  and 
observation  comparisons,  which  are  impossible  to  account  for. 

Figures  3.12  and  3.13  show  the  LWC  comparisons  for 
respectively  the  PS  and  for  the  JW  for  non-convect i ve  (stratus, 
stratocumulus,  and  altocumulus)  clouds.  The  comparison  to  the 
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igure  3.9  j-w  Measured  LWC  vs.  S-F  Predicted  LWC  for 
Stratocumulus  Clouds. 


PARTICLE  SPECTROMETER  MEASURED 


SECTION  4 
CONCLUSIONS 


The  work  and  research  performed  a3  part  of  this  contract  was 
a  logical  extension  of  the  work  and  research  performed  previously 
by  tne  FAA  and  NRL  in  determining  the  character  of  the  atmosphere 
with  regard  to  the  icing  of  aircraft.  The  three  parts  of  the 
work:  collection  of  foreign  data,  determining  the  geographical 

areas  where  icing  would  be  prevalent,  and  finally  to  uevelop  a 
data  collection  scheme  without  aircraft  support  provided  a  logical 
framework  to  develop  a  useful  scheme. 

4.1  DATA  BASE 

The  collection  of  foreign  icing  data  was  nor  as  fruitful  as 
one  would  have  hoped.  There  were  many  disappointments  and 
frustrations.  The  English  in  both  their  meteorological  and 
defense  department  have  data  which  would  have  proved  useful  for 
the  FAA/NRL  data  base.  However,  due  to  events  beyond  our  control, 
the  data  was  not  available  to  us.  However,  we  did  manage  to  ob¬ 
tain  data  from  the  University  of  Manchester  which  was  tower  data 
including  ground  fog.  The  Germans  have  data  from  aircraft  but 
would  only  send  us  a  summary  of  that  data  (3oe  Table  1.2).  The 
Russians  also  have  data  but  no  contact  was  made  with  them.  Those 
arc  the  negative  aspects  of  the  data  collection  effort.  On  the 
pos  .ive  side,  we  obtained  a  large  amount  of  data  from  the  Alborta 
Research  Council  in  Edmonton,  Alberta  and  from  the  Atmospheric 
Environmental  Services  located  in  Ottawa.  Both  of  these  agencies 
perform  cloud  seeding  experiments  on  a  somewhat  steady  plan  and 
have  great  experience  and  expertise  in  the  collection  of  micro- 
physical  data  of  clouds.  Data  obtained  from  the  JTD  Environmental 
Services  Inc.  from  flights  over  the  North  Sea  and  European  coun¬ 
tries  is  also  very  valid  data.  Data  from  the  French  and  South 
Africa  is  precisely  the  same  type  data  as  is  currently  in  the 
FAA/NRL  data  base,  so  that  data  will  be  usoful.  Finally  the  data 
from  Manchester  is  unique  in  that  it  is  tover  data,  including 
ground  fog. 
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Wo  lnoiudad  all  oountrisa  In  our  i ni Mil  request  for  data  bo 
that  each  country  is  aware  of  the  FAA'a  program  in  lolng  data. 

The  Inclusion  of  this  newly  acquired  non-conus  data  Into  the 
praaont  FAA/NRl  will  result  in  a  data  base  that  will  have  data 
principally  from  measurements  made  In  Europe  and  North  America 
along  with  data  from  South  Africa  and  thus  will  be  considerably 
geographically  more  diverse  than  before.  Though  many  or  the 
measurements  ware  made  in  the  Northern  Hemisphere's  severe  lolng 
environment,  Section  2  suggests  that  areas  of  high  frequency  oc¬ 
currence  of  eweru  lolng  were  not  penetrated,  because  these  arese 
have  moderately  heavy  air  traffic,  we  recommend  that  oonal deration 
be  given  for  either  making  additional  flight  measurements  or  im¬ 
plementation  of  the  data  collection  aohsme  discussed  In  Section 
3  . 

Furthei,  son*  >  *  -eo  such  ae  eltoetretua  (Aa),  ai 

tooumulua  (Ao),  on  .uGwQs.nu*  (Ns)  appear  to  be  under- 
represented  in  tnt  now  data  esse.  The  few  modern  measurements 
made  by  the  Air  Force  Geophysloa  Laboratory  show  felrly  respect¬ 
able  LWC  values  (up  to  .4  gm/m*)  for  Ao  end  Ac  so  inclusion  of 
more  as  and  Ao  iceasut  smenta  Into  the  data  base  appears  desirable. 
Ao  end  Ao  asoodletsd  with  strong  wsrrn  fronts  could  present 
horizontally  extensive  regions  of  relatively  high  LWC.  This  type 
of  weather  occurs  often  in  the  mid- 1st  1 tude  westerlies  during  the 
winter. 

Most  of  the  highest  LWC  values  in  the  new  data  base  occur 
within  oonieotlve  clouds  over  the  Sierra  Nevada  range  aru'  over 
northern  Spain.  However,  all  meaeureraento  were  taker.  ,n  winter 
and  spring  months,  arid  since  both  of  these  areas  ars  not  noted  for 
deop  oonvsjt'on  It  la  very  desirable  to  Incorporate  more  warm 
season  deep  convection  mssaur amenta,  An  enoouriflng  flgn  la  the 
recent,  t  riuerporot  Ion  of  Cooperative  Convective  P  reel  p',  tet  Ion 
Experiment  (CCOI’fc)  date,  but  more  meas'irenienv*  such  as  those  ap 
peering  In  the  South  African  end  Canadian  data  ate  needed. 


H,2  THE  ETAC  LHC  DATA  BASE 

The  processing  of  the  ETAC  LWC  (Section  2)  data  base  snows 
that  certain  areaa  In  Asia,  Europe,  and  North  American  are  prone 
to  severe  lolng.  Figure  2,11  shows  that  these  areas  are  not  well 
represented  by  alroraft  measurements  so  any  characterization  of 
the  world-wide  Icing  environment  based  upon  measurements  may  be 
too  lenient.  This  la  especially  true  for  eastern  Asia  and  north¬ 
ern  Europe. 

It  is  obviously  apparent  that  aircraft  measurement  of  all 
areas  suspected  of  having  aavara  icing  la  both  economically  impos¬ 
sible  and  time  consuming.  In  this  modern  era  of  satellite  probes 
and  high  speed  computers  a  better  data  collection  scheme  can  be 
dev  1 aed . 

k.j  ALTERNATE  DATA  COLLECTION  SCHEME 

Based  upon  a  review  of  candidate  modeling  and  remote  sensing 
techniques,  we  conclude  that  the  updated  Siolth-Feddes  model  as 
dlnounsed  In  Section  3  offers  the  boat  hope  for  an  alternative 
dele  development  scheme  to  enhance  the  new  data  base.  As  shown  In 
Sent  Ion  3  tne  updated  Smlth-Feddes  model  gives  s  reliable  predlo- 
tlo  of  LWC  and  temperature  but  poor  predictions  of  the  total  drop 
spectrum.  We  have  proposed  but  not  Implemented  a  method  to  im¬ 
prove  this. 

W n  also  know  that  geographical  location,  in  cloud  position, 
cloud  age  and  interaction  between  olouds  are  important  factors  In¬ 
fluencing  the  drop  size  distribution.  We  believe  that  at  least 
qeogrephy  and  In  oloud  position  could  be  Incorporated  into  the  SFM 
to  yield  drop  elte  dl strl but  Iona .  These  improvements  would  ad¬ 
dress  the  varletlon  of  Cloud  Condensation  Nuclei  (CCN) 
particularly  between  oceanic  and  continental  regions  and  the  shift 
of  the  drop  spectrum  to  higher  radii  with  height  In  clouds.  These 
lmprovementa  could  be  Incorporated  in  a  relatively  simple  manner 
and  validated  agalnet  the  data  base  as  in  Section  3- 
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H.4  CONCLUSIONS 


The  results  of  this  work  have  been  to  add  data  to  the  exist¬ 
ing  data  bank,  to  show  that  some  areas  with  good  potential  for 
hazardous  icing  conditions  have  not  been  subjected  to  aircraft 
measurements ,  and  finally  that  it  is  possible  to  predict  icing 
conditions  to  a  fair  amount  of  reality,  thus  obviating  the  need 
for  costly  and  time  consuming  aircraft  measurements. 
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APPENDIX  A 


This  appendix  lists  all  agencies  to  which  the  inquiry  letter 
concerning  icing  data  was  sent. 

The  response  to  the  inquiry  letter  is  indicated  in  the  response 
line  by  either  (1)  a  "no  response"  meaning  just  that,  (2)  a  "no 
data"  indicating  that  a  letter  was  received  from  the  agency  stating 
that  the  agency  had  no  data  but  with  no  further  information,  or 
(3)  a  letter  was  received  with  information  which  is  summarized  in 
the  response  line. 
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BELGIUM  (Continued) 


04-03 

Association  Beige  des  Pilotes 

et  Navigateurs  de  Ligne  ( ABPNL) 
Avenue  Henri  Dunant  2, 

1140  Bruxelles 

BELGIUM 

(  02)  36  02  64 

FIRST  LETTER  SENT:  8S0712 
RESPONSE:  No  data 


04-04 

Association  Beige  des  Journalistes 
Profess ionnels  de  L ' Aeronaut  ique  et 
de  1' Astronautique  ( ABJPAA) 

Square  de  l'Arbalete  4 
1170  Bruxelles 
BELGIUM 
(02)  73  24  06 

FIRST  LETTER  SENT:  850712 
RESPONSE:  No  data 


04-05 

Mr.  J.  Ooms 
A.B.I  . 

Openlucht  Wandelgang  48 
B  1150  St.  Pxeters-woluwe 
BELGIUM 

FIRST  LETTER  SENT:  850813 
RESPONSE:  No  data 


04-06 

Mr.  DeRidder 
Director  Meter 
Regie  der  Luchtwegen 
Brussels  National  Airport 
B  1930  Zeventa*. 

BELGIUM 

FIRST  LETTER  SENT:  850813 
RESPONSE:  No  data 
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BELGIUM  (Continued) 


0  4-12 

Mr.  Hombla- 
SABENA 

Service  formation  and  Qua!  i  f  icat  ion 
Brussel*  National  Airport 
B1930  Zavantem 
BELGIUM 

011  32  2  720  59  80,  ext.  3966 
FIRST  LETTER  SENT :  850916 

RESPONSE:  No  Response 


04-13 

Mr.  G.  Doumont 

Diractaur  d* Administration, 

R  .  V .  A  * 

Cantra  da  Communication*  Nord, 
Rua  du  Prograt,  84bie,  bta  1 
1000  Bruxalla* 

BELGIUM 

FIRST  LETTER  SENT:  851002 
RESPONSE:  No  rasponaa 


CANADA  (Ganaral ) 

05-01 

Dr.  Larking  Kerwin 

National  Raaaarch  Council  of  Canada 

Montreal  Road 

Ottawa,  Ontario 

CANADA  K1A  0R6 

(  6  13)  993-9  109  (NRC)/(  6  1  3)  9  33-2  37  l  ( St  a  1 1 abra ss ) 
FIRST  LETTER  SENT:  *,850712(cc.  Stallabrass) 
RESPONSE:  No  data 


05-02 

Dr.  Ed  Lozowski 
University  of  Alberta 
Edmonton,  Albarta 
CANADA  T6G  2E8 
( 403)  432-2325  (Univ. ) 

FIRST  LETTER  SENT:  *,850712 

RESPONSE:  Referred  u*  to  Dr.  9.  Humphries  ot  the  Alberta 

Research  Council,  Sea  05-05 
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CANADA  (Continued) 


0  5-07 

Dr.  Robert'  S.  Schemenauer 
Atmospheric  Environment  Service 
4905  Duffer'in  Street 
Downsview,  Ontario 
CANADA  M3H  5T4 
(  4  16)  667-4683 

FIRST  LETTER  SENT:  *,850712 

RESPONSE:  The  AES  sent  data  from  six  projects. 


0  5-08 

Public  Information  Officer 
International  Civil  Aviation 
Organization  (ICAO) 

International  Aviation  Square 
1000  Sherbrooke  St.  W. 

Montreal,  Ouebec 

CANADA  H3A  2R2 

(514)  285-8220 

FIRST  LETTER  SENT:  850712 

RESPONSE:  No  data,  but  recommended  a  Swedish  contact. 


0  5-09 

P.  A.  Corbett,  Exec,  Sec. 
Canadian  Aeronautics  and 
Space  Institute 
#60  75  Sparks  Street 

Ottawa,  Ontario 
CANADA  KIP  5A5 
(613)  234-0191 
FIRST  LETTER  SENT:  850712 
RESPONSE:  No  response 


05-10 

Mr.  Roger  3urgess-Webb 
Managor,  Information  Services 
Canadian  Air  Line  Pilots  Association 
1 300  Steeles  A ve .  E . 

Brampton,  Ontario 
CANADA  L6T  1A2 
(  4  16)  452-8210 
FIRST  LETTER  SENT:  850712 
RESPONSE:  No  response 
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CANADA  (Continued) 

05-15 

Mr.  Ken  Grandia 

Alb.erta  Research  Council 

Atmospheric -Sciences  Department 

7th  Floor,  Terrace  Plaza 

4445  Calgary  Trail  south 

Edmonton,  Alberta 

CANADA  T6H  5R7 

FIRST  LETTER  SENT:  850901 

RESPONSE:  See  05-05 


CANADA  (AES) 

06-01 

Atmospheric  Environmental  Service 
4905  Dufferin  St. 

Downsview,  Ontario 

CANADA  M3H  5T4 

(416)  667-4551 

FIRST  LETTER  SENT:  850712 

RESPONSE:  See  05-06  and  05-07 


06-02 

Dr.  A.  D.  J.  O'Neill 

Director,  Atlantic  Office 

Atmospheric  Environmental  Service 

1496  Bedford  Hwy 

Bedford,  Nova  Scotia 

CANADA  B4A  1E5 

(902)  835-9328 

FIRST  LETTER  SENT:  850712 

RESPONSE:  See  05-06  and  05-07 


06-03 

R.  J.  Fichaud 

Director,  Quebec  Office 

Atmospheric  Env ironmenatl  Service 

100  Alexis  Ntchon  Blvd  -  3rd  Floor 

Ville  St.  Laurent,  Quebec 

CANADA  H4M  2N6 

(514)  333-3000 

FIRST  LETTER  SENT:  850712 

RESPONSE:  Se j  05-06  and  05-07 
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Canada  (Continued) 


06-08 

F.  J.  Lemire 

Canadian  Meteorological  Centre 

2  121  N.  Service  Rd.  #40  4 

Trans  Canada  Hwy 

Dorval,  0C 

CANADA  H9P  1J3 

(  514)  683-7274 

FIRST  LETTER  SENT:  850712 

RESPONSE:  No  reply 


CARIBBEAN 

07-01 

Caribbean  Meteorological  Institute 

P.  0.  Box  130 

Br  idge  town 

BARBADOS 

1-809-425-1362 

FIRST  LETTER  SENT:  850712 

RESPONSE:  No  data 


07-02 

The  Coordinating  Director 

Caribbean  Meteorological  org.  (CMO) 

p.  0.  Box  461 

port-of-Spain 

TRINIDAD- TOBAGO 

1-809-624-3121 

FIRST  LETTER  SENT:  850712 

RESPONSE:  No  data 


07-03 

The  Director 

Meteorological  Services 

Grantley  Adams  International  Airport 

Chri3t  Church 

BARBADOS 

809-428-8521 

FIRST  LETTER  SENT:  850801 
RESPONSE:  No  data 
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FRANCE 


11-01 

Mr.  Martin;  Friedlander 
Centre  d'Aasais  en  Vol  (CEV)  -  France 
French  Military  Testing/Flight  Test  Ctr. 
9  1220  Burt*,  iton  Br i t tany-Fur-Orde , 

FRANCE 

011  33  6084  9570  Ext  3242 
011  33  6171  0249 
FIRST  LETTER  SENT:  *,850712 
RESPONSE:  They  sponsor  Dr.  Gayet 


11-02 

Mr.  J.  F.  (Jean-Franco is )  Gayet 
Laboratoire  Associe  de 
Meteorologie  Physique 
Universite  of  Clermont  II 
B.P.  45,  63170  Aubiere, 

FRANCE 

FIRST  LETTER  SENT:  *,850712 
RESPONSE:  Sent  data  in  the  fall  of  86. 


1 1-03 

Prof.  R.  G.  Soulage 
Universite  de  Clermont  II 
63170  Aubiere 
FRANCE 

FIRST  LETTER  SENT:  *,850712 

RESPONSE:  He  is  Dr.  Gayet's  supervisor. 


1 1-04 

Madam  H.  Bouilloud 
Ministry  of  Defense 
STPA/CIN 

4  Avenue  Deiaporte  d'essy 
75996  Paris  Armees 
FRANCE 
1  552  5319 

FIRST  LETTER  SENT:  850712 
RESPONSE:  They  sponsor  Dr.  Gayet 
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FRANCE  (Continued) 

11-10 

Meteorologie  Nationals 
Respondeur  Renseignements  lie  de 
France  et -Normandie 
2  av  Rapp  7 
Paris , 

FRANCE 

(1)  S55.95.90 

FIRST  LETTER  SENT:  350712 
RESPONSE:  No  response 


11-11 

Mission  des  etudes  et  de  la  recherche 
14,  boulevard  du  General  Leclerc 
92524  Neui  l  ley-sur-Se  me  Cedex 
Paris , 

FRANCE 

(1)  756.12.12 

FIRST  LETTER  SENT:  050712 
RESPONSE:  No  response 


11-12 

Centre  national  de  la 
recherche  sc  lent  ?.  f  lque 
15,  quai  Anatale  France  7 
Paris , 

FRANCE 

(1)  555.92.25 

FIRST  LETTER  SENT:  850712 
RESPONSE:  No  response 


GERMANY 

12-01 

Dr.  Hans  Eberhard  Hoffmann 
Institute  of  Atmospheric  Physics 
( Institut  fur  Physik  der  Atrosphare) 

Deutsche  Forschungs-  und  Ver suchansta 1 t 
fur  Luft-  und  Raumfahrt  (DFVLR) 

Oberpfaf fenhofen  D-8031 
8031  We ss ling 
GERMANY  FR 

011  4981  5328  Ext  579  (081  53/28-579) 

FIRST  LETTER  SENT:  850712 

RESPONSE:  Dr.  Hoffmann  has  sent  us  summaries  of  his  data.  For 

reasons  not  clear  to  UDRI,  he  did  not  send  UDRI  the  magnetic 
tape  data. 
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Germany  (Continued) 


12-06 

Mr.  peter  Schramm 
German  Defense  Test  Center 
or  German  Experimental  Center 
( Bondeswehr  Erprobungsstelle ) 

Er probung ss telle ,  61  der  Bundeswehr 
Dezernat  234  Flugplatz, 

8072  Manching 
WEST  GERMANY 
08459-802174 

FIRST  LETTER  3ENT :  850712 

RESPONSE:  See  12-01 


12-07 

Mr.  Kurt  Uwira 
German  Defense  Test  Center 
or  German  Experimental  Center 
(Bondeswehr  Erprobungsstelle) 
Bundesamt  fur  Wehrtechnik  und 
Beschaffung  AFB  LG  III 
Flugplatz 
8072  Manching 
WEST  GERMANY 

FIRST  LETTER  SENT:  850712 
RESPONSE:  See  12-01 


12-08 

Herr  Richter 

German  Defense  Test  Center 
or  German  Experimental  Center 
Bundesamt  fur  wehrtechnik  und  Beshaffmg 
Erprobungsstelle,  61  der  Bundeswehr 
8072  Manching 
WEST  GERMANY 

FIRST  LETTER  SENT:  850712(copy  of  12-06) 
RESPONSE:  See  12-01 


12-09 

Mr.  Eckert 

German  Defense  Test  Center 
or  German  Experimental  Center 
Bundeminister  fur  Verteidigung 
RU  IV  6,  F.  Hel  Hern  Eckert 
53  Bonn 
WEST  GERMANY 

FIRST  LETTER  SENT:  8S0712(copy  of  12-06) 
RESPONSE:  See  12-01 
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ITALY 


15-01 

Associazione  Industrie  Aerospazial  i 
via  Naziovale  200 
Roma , 

ITALY 

(06)  46  02  47 

FIRST  LETTER  SENT:  850712 
RESPONSE:  No  response 


15-02 

AO PA  Italia 

Organ i zzoz ione  dell  'Aviazione 
private  e  d'Affan 
Corso  Magenta  56, 

20123  Milano, 

ITALY 

(  02)  87  38  02 

FIRST  LETTER  SENT:  8S0712 
RESPONSE:  No  response 


15-03 

Associazione  Nazionale  Piloti 
Aviazione  Commerc lale  (ANPAC) 
Viale  dell ' Esperanto  71 
00144  Roma 
ITALY 

(06)  591  04  11 

FIRST  LETTER  SENT:  850712 

RESPONSE:  No  response 


15-04 

Associazione  Nazionale  Piloti 
Aviazione  Generale  ( AN  PAG ) 
Via  Mam i am  15/2 
20127  Milano 
ITALY 

(02)  289  60  59 

FIRST  LETTER  SENT:  850712 

RESPONSE:  No  response 


15-05 

Unione  Giornalisti  Aeronautici  Italim 
Corso  Trieste  10 
00198  Roma 
ITALY 

(06)  85  51  71 

FIRST  LETTER  SENT:  850712 
RESPONSE:  No  response 
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THE  NETHERLANDS 


19-02 

prins  Maurit3  Lab.,  TNO 

(Dutch  org.  foe  Applied  Research) 

Lande  Kleiweg  137 

228  GJ  Rijswijk 

THE  NETHERLANDS 

FIRST  LETTER  SENT:  850712 

RESPONSE:  Sent  letter  to  the  Dutch  National  Air  and  Space 

Laboratory  who  in  turn  sent  it  to  Fokker  Aircraft.  See  19-08 
of  this  Appendix. 


19-03 

AOPA  Netherlands 

jozef  israelsplein  8 

De  n  Ha  ag 

THE  NETHERLANDS 

FIRST  LETTER  SENT:  850712 

RESPONSE:  No  response 


19-04 

Vereniging  van  Nederlandse 
Verkeersvl iegers  (VNV) 

Dutch  Airline  Pilots  Association 
Ams  terdamseweg  138 
Amstelveen , 

THE  NETHERLANDS 
(020)  41  05  55 
FIRST  LETTER  SENT:  850712 
RESPONSE:  No  response 


19-05 

Nederlandse  Vereniging  van  Lucht-  en 
Ruimtevaart-  Publicisten  (NVLP) 
Jozef  israelsplein  8 
Den  Haag 
THE  NETHERLANDS 
(070)  24  72  52 
FIRST  LETTER  SENT:  850712 
RESPONSE:  No  response 


19-06 

Dr.  E.  Hosstee 

Meteo  Amsterdam  Airport 

THE  NETHERLANDS 

FIRST  LETTER  SENT:  850712 

RESPONSE?  No  response 
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PORTUGAL 


22-01 

Sindicato  Nacional  de  pessoal 
da  Voo  d*  Aviacao  Civil  (SNPVAC' 
pcaca  Pasteur  11 
R/C-D,  Li3boa 
PORTUGAL 
72  87  74 

FIRST  LETTER  SENT:  850712 
RESPONSE:  No  response 


SCOTLAND 

23-01 

Scottish  Meteorological  Society 
Ed  inburgh, 

SCOTLAND 

FIRST  LETTER  SENT:  850712 
RESPONSE:  No  response 


SOUTH  AFRICA 

24-01 

Roelof  Bruintjes 
South  African  weather  Bureau 
P.  O.  Box  169 
Irene,  1675 

REPUBLIC  OF  SOUTH  AFRICA 
011-27-12-290-000  (6  hours  ahead  of  U.S.) 

FIRST  LETTER  SENT:  *,850712 

RESPONSE:  Received  two  magnetic  tapes  of  data. 


24-02 

Council  for  Scientific  and  Industrial 
Research  (CSIR) 

CSIR/WNNR 

Attn:  Cloud  Physics  Group 
Pretoria , 

REPUBLIC  OF  SOUTH  AFRICA 
FIRST  LETTER  SENT:  850712 

RESPONSE:  Suggested  24-03  ,  24-0  4 ,  and  24-05 
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SPAIN  (Continued) 


25-02 

Ag rupee ion  Sindical  de 

Pilotes  de  Lineas  Aereas  (ASPLA) 
paaeo  del  Pfado  18-20 
Planta  5a 
Madrid  14 
SPAIN 

FIRST  LETTER  SENT:  850712 
RESPONSE:  No  response 


SWEDEN 

26-02 

AO PA  Sweden 

Svenska  Allmanf lygforemng 
Fack,  161  10  Bromma  10 
SWEDEN 

(08)  29  50  00 

FIRST  LETTER  SENT:  850712 

RESPONSE:  No  data.  They  also  indicated  that  neither  the 

Swedish  Board  of  Civil  Aviation,  the  Swedish  Meteorological  and 
Hydrological  institute,  the  Stockholm  University,  the  Aero¬ 
nautical  Research  institute  of  Sweden,  nor  the  Royal  Technical 
High  School  of  Stockholm  had  any  data. 


26-03 

Svensk  Pilotforening  (SPF) 

Swedish  Airline  Pilots  Association 
Olofsgatan  10,  2  tr, 

111  36  Stockholm  C. 

SWEDEN 

(08)  10  94  34 

FIRST  LETTER  SENT:  850712 
RESPONSE:  No  response 


26-04 

John  Ogren 

University  of  Stockholm 

Dept,  of  Meteorology/ Arrhenius  Lab. 

S- 106  91, 

Stockholm, 

SWEDEN 

011-46-8-162000  (Univ.  of  Stockholm) 
FIRST  LETTER  SENT:  850712 
RESPONSE:  No  data 

Mr.  Heintzenberg 
RESPONSE:  No  response 
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SWITZERLAND  (Continued) 


27-02 

Global  Atmospheric  Research 
Programme  (GARP) 

GARP  Activities  office 
c/o  WHO 

Case  Postale  5,  CH  1211 
Geneve  20, 

SWITZERLAND 
34  64  00 

FIRST  LETTER  SENT:  850712 
RESPONSE:  No  response 


27-03 

Vereinigung  der  Schweizereschen 
Fluggeug  industr ie 
Association  Suisse  de  1* Industrie 
Aeronautique  (ASIA)  Theaterplatz  4, 
5400  Baden, 

SWITZERLAND 
(056)  2  30  90 

FIRST  LETTER  SENT:  850712 
RESPONSE:  No  response 


27-04 

AOPA  Switzerland 
P.  O.  Box  151 
8058  Zurich  Airport 
Zurich, 

SWITZERLAND 
(051)  84  01  85 
FIRST  LETTER  SENT:  850712 
RESPONSE:  No  data 


27-05 

Aeropers  (Schweiz.) 

8152  Glattbrugg, 
Rietstrasse  17 
SWITZERLAND 

FIRST  LETTER  SENT:  850712 
RESPONSE:  No  response 


27-06 

Schwe izer ische  Meteorologische  Anstalt 
Swiss  Meteorological  Institution 
Krahbuhlstrasse  58 
CH-8044  Zurich  SWITZERLAND 
FIRST  LETTER  SENT:  850807 
RESPONSE:  No  response 
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UNITED  KINGDOM  (Continued) 


29-04 

Dr.  J.  T.  Cansdale 

Royal  Aircraft  Establishment 

Farnsborough ,  Harts., 

ENGLAND 

011  44  252  24461  Ext.  2491 

FIRST  LETTER  SENT:  *,350712 

RESPONSE:  Sent  letter  to  A&AEE  (29-01) 


29-05 

G.  M .  E.  white 
Secretary 

Civil  Aviation  Authority 
CAA  House,  45-59  Kingsway 
London 

UNITED  KINGDOM  WC2B  6TE 
01-379  7311 

FIRST  LETTER  SENT:  350712 
RESPONSE:  No  data 


29-07 

Frank  Atkinson 
British  Airways  Helicopters 
Gatwick  Airport 
Surrey , 

ENGLAND 

FIRST  LETTER  SENT:  850712 
RESPONSE:  No  response 


29-09 

Peter  Soliz 
Major  US AF 

Chief,  Geophysics  and  Space 
European  Office  of  Aerospace  R&D 
2  23/231  Old  Marylebone  Rd . 
London,  NW1  5th 
ENGLAND 

(01)  409-4437  TELEX-299739 

FIRST  LETTER  SENT:  850712 
RESPONSE:  See  29-01 
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UNITED  KINGDOM  (Continued) 


29-15 

Mr.  R.  0.  Belton,  Secretary 
Board  Air  Safety  Review  Committee 
British  Airways, 

Heathrow  Airport 
P.  0.  Box  10 
Hounslow , 

ENGLAND  TW6  2JA 
011-44-01-759-5511 
FIRST  LETTER  SENT :  8507  12 

RESPONSE:  No  data 


29-16 

Mr.  Alan  A.  Woodfield, 

Royal  Aircraft  Establishment 
General  Aerodynamics  Section, 
Flight  Research  Division 
Bedford , 

ENGLAND 

011-44-234-55241 

FIRST  LETTER  SENT:  850712 

RESPONSE:  No  response 


29-17 

Capt.  R.  D.  Hillary 
Commonwealth  Advisory  Aeronautical 
Research  Council  (CAARC) 
c/o  Ministry  of  Defense 
Old  War  Office  Bldg.  -  Room  SW1A 
Whitehall,  London, 

UNITED  KINGDOM  SW1A  2EU 
218  0838 

FIRST  LETTER  SENT:  850712 
RESPONSE :  No  data 


29-18 

Rolls-Royce  Ltd. 

65  Buckingham  Gate 
London,  SW1, 

UNITED  KINGDOM 

FIRST  LETTER  SENT:  850712 

RESPONSE:  No  response 
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UNITED  KINGDOM  (Continued) 

Agency: 

The  University  of  Manchester 
Institute  of  Science  and  Technology 

Contact :  • 

Professor  John  Latham 
The  University  of  Manchester 
Institute  of  Science  and  Technology 
P.  0.  Box  88 
Manchester  M60  19D 

061-236-3311 

RESPONSE:  Received  five  days  of  data. 


USSR 

30-01 

prof.  I.  P*  Mann 

Soviet  Geophysical  Committee 

Molodezhnaya  3, 

Moscow  117296, 

USSR 

Not  available 

FIRST  LETTER  SENT:  850712 
RESPONSE:  No  response 


30-02 

The  Institute  of  Atmospheric  Physics 
University  of  Moscow 
Leninsky  Gory 
Moscow , 

USSR 

FIRST  LETTER  SENT:  850712 
RESPONSE:  No  response 


30-03 

Dr.  Leonid  T.  Matveev 

Leningrad  Hydrometeorological  Institute 

Leningrad 

USSR 

FIRST  LETTER  SENT:  831105 
RESPONSE:  No  response 
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APPENDIX  B 


The  following  pages  comprise  the  mailing  that  was  sent  to  each 
of  the  agencies  listed  in  Appendix  A. 
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The  University  of  Dayton 


( Date ) 


(Contact's  Name) 
(Agency  Name) 
(Agency  Address) 


Dear  (Contact's  name): 

The  Federal  Aviation  Administration  (FAA)  recently 
completed  an  aircraft  icing  data  gathering  and  analysis  effort 
which  resulted  in  a  new  characterization  of  the  aircraft  icing 
environment  in  supercooled  clouds  up  to  3  tan  above  ground  level 
for  the  continental  United  States.  The  FAA  is  now  continuing  to 
characterize  the  aircraft  icing  atmosphere  for  all  conditions  and 
for  all  altitudes  both  in  the  U.S.A.  and  worldwide  in  a  planned 
sequence . 

The  Icing  data  base  is  being  developed  and  maintained  at 
the  U.S.  Naval  Research  Laboratory  under  tne  direction  of  Dr. 
Richard  Jeck  of  NRL  and  sponsored  by  the  FAA. 

The  role  of  UDRI  as  explained  in  the  cover  letter  is  to 
establish  contact  with  governments  and  scientific  groups  outside 
the  U.S.A.  and  arrange  for  the  retrieval  of  pertinent  aircraft 
icing  data  so  that  when  added  to  the  already  existing  continental 
U.S.A.  icing  data  base,  a  truly  world-wide  data  base  will  be 
available  to  all  participating  countries. 

We  are  interested  in  retrieving  the  following  types  of  data 
which  we  have  labeled  Primary,  Secondary,  and  Supplementary. 

Primary  Data  (Required  to  characterize  the  icing  environment) 

Simultaneous.  measurements  in  ground  fog  or  in/beneath 
clouds  containing  any  supercooled  water  droplets  including 
freezing  precipitation  or  ice  crystals  including  snow,  or  com¬ 
binations  thereof,  of  the  following  variables: 
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RESEARCH  INSTITUTE 
Applied  Syttemt  Analyst 
300  Col  lest  Park  Dayton,  Ohio  43469-0001 


Variable 


Liquid  Water  Content  (or  -  redundant  measurements 

equivalent)  preferred 

Droplet  Median  volume  Diameter 


Air  Temperature  -  corrected  preferred 

Aircraft  Altitude  -  pressure  altitude-  and 

above  ground  level 
preferred 


Aircraft  True  Airspeed 

Time/Duration/Distance  of 
Icing  Event 

Secondary  Data 

Secondary  data  is  desirable  but  not  necessary, 
Variable 


Droplet  Size  Spectra 
Droplet  Concentration 
Aircraft  Location 
Aircraft  Attitude 

Terrain  Height 
Supplementary  Information 

Supplementary  information 


-  distribution 

-  number  density 
•  geographic 

-  climb,  level  flight, 
descent,  etc. 

-  mean  sea  level  elevation 

s  highly  desirable. 


Technical  reports  on  data  analysis  or  case  studies. 

Summary  reports  on  flight  operation  or  field  projects. 

Scientific  observation  notes  or  flight  logs  from  pilot  or 
onboard  scientist,  especially  information  on  cloud 
conditions,  icing  conditions,  precipitation  type  and 
intensity,  instrument  performance. 

Information  on  data  quality,  error  correction  procedures. 


Note:  Preferred  data  medium  is  digital  computer  tapes  with 
documented  format. 
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With  this  information  in  mind  we  would  like  to  ask  you  the 
following > 

1.  Do  any  of  the  above  mentioned  kinds  of  data  exist  at  your 
agency? 

2.  If  so,  19  it  possible  for  L'DRI  to  obtain  the  data  for  the 
FAA? 

3.  If  the  data  does  exist,  in  what  form  does  it  exist,  i.e., 
magnetic  tape,  reports,  computer  printout,  etc,  (Please  be  as 
specific  as  possible)  and  is  it  possible  to  send  samples  of 
the  data? 

4.  Would  it  be  possible  for  you  to  send  us  a  list  of  airborne 
cloud  physics  field  projects  over  the  past  ten  years  end  a 
lUt  of  references  to  related  technical  reports  (published  or 
internal),  conference  papers,  and  published  scientific 
articles,  that  your  agency  has  been  involved  in? 

5.  Could  you  recommend  other  places  to  inquire  among 
universities,  military,  and  government  agencies,  either  in 
your  country  or  otners? 

6*  Could  you  furnish  us  with  the  name(s)  of  the  proper 
person(s)  to  contact  on  this  matter  and  their  correct 
address  and  phone  number? 

The  UDRI  and  the  FAA  both  recognize  and  appreciate  any 
efforts  you  may  be  able  to  extend  to  us.  The  FAA  has  already 
stated  that  they  consider  the  creation  of  this  world-wide  data 
base  to  be  an  international  effort  and  as  su  :h  are  willing  to 
share  the  completed  data  base  with  all  participants. 

Thank  you  for  your  time  and  consideration  in  this  important 
endeavor , 

With  best  regards*, 

Nicholas  A,  (Tf^Jl'sr 

Senior  Research  physicist 

NAE/gv 
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A"*«*  Crv 

N*'"  Jtfity  08*05 


To  whom  it  may  concerns 

The  University  of  Dayton  is  currently  under  contract  with 
the  Federal  Aviation  Administration  (FAA) ,  to  conduct  a 
research  effort  to  better  define  the  atmospheric  parameters 
conducive  to  the  formation  of  in-flight  icing  on  aircraft. 

The  results  of  the  effort  will  be  usee1  to  aid  in  the  design 
of  aircraft  anti-icing  and  deicing  equipment. 

The  University  of  Dayton,  intends  to  contact  various  foreign 
cloud  physics  groups  and  other  related  activities,  concerning 
information  on  this  matter.  Your  cooperation  with  the 
University  of  Dayton  is  greatly  appreciated. 

Sincerely, 

FAA  Contracting  Officer 


V  ~4S—" 

Ddnna  C.  Land 
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Measured  LWC  vs.  S-F  Predicted  LWC  for  Cumulus 
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Figure  3.12  P-S  Measured  LWC  vs.  S-F  Predicted  LWC  for  ST,  SC,  and 
AC  Cloud  Types. 
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particle  spectrometer  is  particularly  good,  we  note  that  the  com¬ 
parison  continues  to  be  good  even  for  higher  LWC  values 
(admittedly  too  few). 

Finally  Figures  3.14  and  3.15  show  the  LWC  com¬ 
parisons  between  the  JW  and  the  PS  for  respectively  non-convective 
and  convective  clouds.  The  percent  error  for  non-convective 
clouds  is  greater  than  that  for  the  comparison  between  the  PS  and 
the  Smith-Feddes  model.  It  is  interesting  that  the  LWC  measure¬ 
ments  for  cumulus  have  a  smaller  percent  error  between  the  two 
devices  than  between  the  model  and  the  LWC  measurement  devices. 

It  is  tempting  to  speculate  that  the  percent  error  is  relatively 
larger  for  smaller  LWC  values  and  is  expected  to  decrease  with  in¬ 
creasing  LWC.  Some  of  the  discrepancy  between  the  model  and 
measurements  can  be  attributed  to  the  sampling  strategy. 

3.5  CONCLUSIONS  FROM  THE  VALIDATION  STUDY 

(1)  The  model  compares  slightly  better  to  the  particle 
spectrometer  than  to  the  Johnson-Williams  probe. 

(2)  The  model's  percent  error  is  about  the  same  for  all 
types  of  clouds,  perhaps  slightly  better  for  Cu  than 
for  other  types  of  clouds. 

(3)  The  model's  predictive  error  is  no  worse  than  the 
measuring  devices'  error. 

( *0  Even  with  an  approximate  -10?  error,  the  uncertainty  in 
the  icing  severity  is  minimal.  For  example,  given  a 
moderate  to  severe  icing  condition,  which  has  a  mean 
volume  drop  diameter  of  15  microns,  and  a  LWC  of  1 

3 

gtn/m  ,  the  40$  uncertainty  would  produce  a  LWC  range  of 
0.6  (moderate)  to  1.4  gm/m^  (severe)  icing  condition. 
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gure  3.14  J-W  Measured  LWC  vs.  P-S  Measured  LWC  for  ST,  SC,  and 
AC  Cloud  Types. 
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3.6  APPLICATION  OF  THE  SMITH-FEDDES  MODEL  TO  ALTERNATE  DATA 

DEVELOPMENT 

With  improvements  to  the  drop  size  calculation  as  mentioned 
above,  the  modified  updated  Sraith-Feddes  model  satisfies  the 
criteria  for  an  alternate  data  development  scheme  i.e.,  it  gives 
useful  values  of  at  least  LWC  and  temperature  and  reasonable 
values  of  MVD.  In  fact,  the  agreement  of  JW  and  PS  measurements 
is  not  appreciably  better  than  for  the  measurements  and  the  model 
calculations.  With  such  a  tool  available,  the  question  becomes 
one  of  use;  how  should  this  version  of  the  Smith-Feddes  model  be 
used  to  generate  additional  data  to  extend  the  validity  of  the 
FAA/NRL  data  base?  We  envision  two  possible  strategies  to  do 
this.  One  strategy  is  to  make  a  statistically  significant  number 
of  Smith-Feddes  calculations  for  data  sparse  regions  where  the 
icing  environment  has  been  identified  as  severe  (Section  2).  The 
second  strategy  is  to  make  a  parametric  study  by  varying  the 
model’s  input  and  then  to  assess  the  climatological  likelihood  of 
input  conditions  that  predict  severe  icing. 

Under  the  first  strategy,  the  Smith-Feddes  model  would  be 
run  for  a  large  number  of  data  inputs.  The  data  would  arise  from 
RTNEPH  and  AFGWC  temperature  and  height  analyses  for  the 
geographical  regions  identified  in  Section  2  as  prone  to  severe 
icing.  A  major  question  is  the  quality  of  the  RTNEPH  analysis  for 
the  regions  identified  in  Section  2  as  well  as  the  vertical 
resolution  of  temperature  and  height  analyses.  A  final  question 
is  how  many. runs  are  enough?  Perhaps  the  guidelines  set  forth  by 
Jeck  (1983)  for  data  miles  for  various  weather  conditions  can  be 
adapted  to  the  alternate  data  application.  Otherwise  we  could  use 
one  year's  data  of  twice  per  day  observations. 

For  the  second  strategy,  a  parametric  variation  of  Smith- 
Feddes  input  over  conceivable  ranges  would  give  a  number  of  output 
scenarios.  It  would  only  be  necessary  to  consider  those  scenarios 
that  predict  severe  icing.  We  would  then  establish  the 
climatological  llklihood  of  the  severe  icing  scenarios  for  the 
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geographical  regions  identified  in  Section  2.  This  would  estab¬ 
lish  the  probability  of  severe  icing  for  these  regions.  By 
performing  the  identical  procedure  for  areas  of  severe  icing  for 
which  measurements  exist  we  can  establish  the  validity  of  this  ap¬ 
proach  and  establish  whether  or  not  the  data  base  is 
representative  of  the  world-wide  icing  environment. 
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